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ABSTRACT 

We present numerical simulations of galaxy clusters with stochastic heating from active galac- 
tic nuclei (AGN) that are able to reproduce the observed entropy and temperature profiles of 
non-cool-core (NCC) clusters. Our study uses TV -body hydrodynamical simulations to inves- 
tigate how star formation, metal production, black hole accretion, and the associated feedback 
from supernovae and AGN, heat and enrich diffuse gas in galaxy clusters. We assess how 
different implementations of these processes affect the thermal and chemical properties of the 
intracluster medium (ICM), using high-quality X-ray observations of local clusters to con- 
strain our models. For the purposes of this study we have resimulated a sample of 25 massive 
galaxy clusters extracted from the Millennium Simulation. Sub-grid physics is handled using 
a semi-analytic model of galaxy formation, thus guaranteeing that the source of feedback in 
our simulations is a population of galaxies with realistic properties. We find that supernova 
feedback has no effect on the entropy and metallicity structure of the ICM, regardless of the 
method used to inject energy and metals into the diffuse gas. By including AGN feedback, 
we are able to explain the observed entropy and metallicity profiles of clusters, as well as the 
X-ray luminosity-temperature scaling relation for NCC systems. A stochastic model of AGN 
energy injection motivated by anisotropic jet heating - presented for the first time here - is 
crucial for this success. 

With the addition of metal-dependent radiative cooling, our model is also able to produce 
CC clusters, without over-cooling of gas in dense, central regions. 

Key words: hydrodynamics - methods: N-body simulations - galaxies: clusters: general - 
galaxies: cooling flows - X-rays: galaxies: clusters. 



1 INTRODUCTION 
1.1 Background 

Clusters of galaxies are believed to be the largest gravitationally- 
bound objects in the Universe. Their deep gravitational potential 
well means that the largest clusters are 'closed boxes', in the sense 
that baryons ejected from cluster galaxies by supernova (SN) ex- 
plosions and active galactic nuclei (AGN) do not escape the clus- 
ter completely, but instead end up in the hot, diffuse plasma that 
fills the space between cluster galaxies - the intracluster medium 
(ICM). 

The thermal properties of intracluster gas, which can be mea- 
sured with X-ray telescopes such as Chandra, XMM-Newton and 
Suzaku, thus provide a unique fossil record of the physical pro- 
cesses important in galaxy and galaxy cluster formation and evolu- 
tion, such as radiative cooling, star formation, black hole accretion 
and the subsequent feedback from supernovae (SNe) and AGN. In 
addition, measurements of the ICM chemical abundances yield in- 
formation about the production of heavy elements in stars in mem- 
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ber galaxies, providing constraints on nucleosynthesis, and the pro- 
cesses responsible for their transport into the ICM. 

A key diagnostic of the thermal state of intracluster gas is pro- 
vided by the gas entrop}Q. Entropy remains unchanged under adi- 
abatic processes, such as gravitational compression, but increases 
when heat energy is introduced and decreases when radiative cool- 
ing carries heat energy away, thus providing an indicator of the 
non-gravitational processes important in cluster formation. 

In recent years, spatially-resolved observations have facili- 
tated a detailed examination of the radial distribution of entropy 
in clusters. Observed entropy profiles are typically found to scale 



as if a 



at large cluster-centric radii, r > 0.1r2 U (fl 



(e.g. IPonman et al.l 120031: Isun et alj|2009l ; ICavagnol o et al. 1 120091 : 
Sander soiTet alj|2009l ; IPratt et alj|20ld) . This power-law scaling 
agrees with that predicted by simple analytical models based on 



1 We define the gas entropy as K = k^T/n2~ , where kg is Boltz- 
mann's constant, T is the gas temperature, n c is the electron number density 
and 7 = 5/3 is the ratio of specific heats for a monoatomic ideal gas. 

2 We define r a as the radius of a spherical volume within which the mean 
matter density is A times the critical density at the redshift of interest. The 
mass enclosed within this sphere is denoted by Ma ■ 
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spherical collapse jTozzi & Normanll200lh and cosmological hy- 
drodynamical simulations that include gravitational shock heat- 
ing only (e.g. lVoit et alJl2005l : lNagai et alj|2007l : lshort et alJl2010L 
hereafter STY 10). 

However, in the inner regions of clusters, observations have 
unveiled the presence of a mass-dependent entropy excess with re- 
spect to theoretical expectations, and a large dispersion in central 
entropy values (e.g. IP onman et al. 1999; Llovd-Davies et al. 2000; 



Ponm an et al.l 120031: IPratt et alj [2006 ; iMorandi & Ettoril \2001 : 



Cavagnolo et al .l2009l ; |Prattetal.l201()h . The source of this entropy 



excess is likely to be a combination of non-gravitational heating 
from astrophysical sources, such as SNe and AGN, and cooling 
processes. There is evidence that the distribution of central entropy 
values is bimodal, with morphologically disturbed non-cool-core 
(NCC) systems having an elevated central entropy compared to dy- 
namic ally relaxed cool-core (CC) systems (e.g. ICavagnolo et al.l 
2009). It is thought that the association of unrelaxed morphology 
with a high central entropy is an indication that either cool cores 
are destroyed by mergers, or that cool cores have never been able 
to form in these objects. 

The chemical state of the ICM is characterised by its metallic- 
ity, the proportion of chemical elements present heavier than H and 
He. X-ray spectroscopy of galaxy clusters has revealed emission 
features from a variety of chemical elements, O, Ne, Mg, Si, S, Ar, 
Ca, Fe and Ni, all of which are synthesised in stars and transported 
to the ICM by processes such as ram-pressure stripping, galactic 
winds and AGN outflows. Type la SNe produce a large amount of 
Fe, Ni, Si, S, Ar and Ca, but compared to Type II SNe, they only 
produce very small amounts of O, Ne and Mg. Type la SN products 
are found to dominate in clus ter cores, whereas Type II SN products 
are m ore evenly distributed jFinoguenov et al. 2000; Tamura et al.l 
l200ll : seejBohringer & Wernedl2010] for a recent review). This can 
be explained by early homogeneous enrichment by Type II SNe, 
which produce a-elements in the protocluster phase, and a subse- 
quent, more centrally-peaked enrichment by Type la SNe, which 
have longer delay times and continue to explode in the central cD 
galaxy long after the cluster is formed. 

The first detailed measurements of spatia l abun dance distri- 
butions were made by |Pe G randi & Molend| ( l200ll) . using data 
from BeppoSAX, who measured the radial Fe abundance pro- 
files for a sample of massive clusters. They found that CC clus- 
ters have a sharp Fe abundance peak in central regions, whereas 
NCC clusters have flat Fe abundance profiles. Subsequent obser- 
vations with XMM-Newton and Chandra have confirmed this di- 
chot omy between the me t allicity distribution in CC and NCC clus- 
ters jTamura et alj |2004 IVikhlinin et al.l 120051: IPratt et alj 12007 



Baldi et al. 2007; lMaughan et alJl2008l:lLeccardi & Molendil 



2008 



Matsushitall201ll) . 

It is thought that cluster mergers, and the subsequent mixing 
of intracluster gas, are responsible for destroying the central abun- 
dance peak found in CC clusters. However, some systems with a 
highly disturbed morphology are als o found to have a high cen- 
tral m etallicity. iLeccardi et alj d201 61 see also|Rossetti & Molendil 
l20ld) suggest that these objects correspond to relaxed CC systems 
that have undergone a major merger, or a significant AGN heating 
event, very recently, so that mixing processes have not yet had suffi- 
cient time to fully erase low-entropy gas and the central abundance 
peak. 

Explaining the observed thermal and chemical properties of 
the ICM from a theoretical perspective requires a detailed un- 
derstanding of the complex interplay between large-scale gravita- 
tional dynamics and the various small-scale astrophysical processes 



mentioned above. Numerical cosmological hydrodynamical simu- 
lations have emerged as the primary tool with which to tackle this 
problem. There has been considerable effort to include the pro- 
cesses relevant for cluster forma tion and evolution in simul ations 
in a self-consistent manner; see iBorgani & Kravtsovl (2009) for a 
recent review. However, an explicit treatment is unfeasible since 
these processes all occur on scales much smaller than can be re- 
solved with present computational resources. 

Hydrodynamical simulations that include models of radiative 
cooling, star formation, metal production and galactic winds gen- 
erally fail to reproduce observed ICM temperature, entropy and 
metallicity profiles. Simulated temperature profiles typically have 
a sharp spike at small cluster-centric radii, followed by a rapid 
drop i n temperature moving f urther into the core ( e. g. IVald arnini 
20031 iTornatore et alj 120031: IBorgani et alj 12004 ; borneo et all 
2006l : ISiiacki et alj|2007l ; iNagai et alj|2007h , in clear conflict with 



the smo othly declining (flat) profiles of observed CC (NCC) clus- 
ters (e .g. lSanderson et alj2006l : rvikhlinin et alj2006l ; lArnaud et all 
2010). This is due to the adiabatic compression of gas flowing 
in from cluster outskirts to maintain pressure support, following 
too much gas cooling out of the hot phase. This over-cooling 
causes excessive star formation in cluster cores, with predicted stel- 
lar fractions bei ng about a factor of 2 larger than the observed 
value of ~ 10% dBalogh et al.ll200ll ; lLin et alJl2003l ; lBalogh et alj 
2008), which in turn leads to excessive Fe production in cen- 
tral re gions, generatin g steeper abundance profiles than observed 



(e.g. IValdaminjll2003|; ITornatore et alj|2004l; iRomeo et alj 120061: 



ITornatore et al.ll2007l ; lDave et alj|2008h 

It is generally accepted that the solution to the over-cooling 
problem in hydrodynamical simulations is extra heat input from 
AGN. Simple analytical arguments convincingly show that the en- 
ergy liberated by accretion onto a central super-massive black hole 
is sufficient to suppress gas cooling and thus quench star formation. 
The precise details of how this energy is transferred to the ICM are 
not well understood at present, but it appears that there are two ma- 
jor channels via which black hol es interact with their surroundings 
(see lMcNamara &~N ulsen 2007 for a review). 

At high redshift, mergers of gas-rich galaxies occur fre- 
quently and are expected to funnel copious amounts of cold gas 
towards galactic centres, leading to high black hole accretion rates 
and radiating enough energy to support the luminosities of pow- 
erful quasars. Quasar-induced outflows hav e been observation- 
ally confirmed in a number of ca ses (e. g. IChartas et alj 120031: 
[ Crenshaw et alj|2003l: IPounds et all 120031 : Banguly & Brothertonl 
l2008l : lDunnetal.ll2010l) . 

Evidence for another mode of AGN feedback, not related 
to quasar activity, can be seen in nearby CC clusters, which of- 
ten contain radio-loud X-ray cavities in the ICM. These bubbles 
are thought to be inflated by re lativistic jets launched from the 
centra l super-massive bl ack hole (Bla nton et ai]|200ll ; Birzan et al 



20041; iMcNamara et aU 120051; iFabian et alj 120061; iMorita et al 
20061: IJetha et alJ|2008l; iGastaldello et alj|2009l ; iDong et alJl2oioF 
Giacintucci et alj I2OIII) . Bubbles may rise buoyantly, remov- 



ing some of the central cool, enriched gas and allowing it to 
mix with hotter gas in the outer regions of groups and clus- 
ters. Together with the accompanying mechanical heating, this 
can constitute an efficient mechanism for suppressing cool- 
ing flows, and redistributing metals throughout the ICM. Such 
flows are seen in simulations of idealis ed clusters, performed 



with hydrodynamical mesh codes (e.g. [ churazov et alj [200 ll: 
Ouilis et al .1 [200ll : iRuszkowski & Begelmanll2002l ; iBruggen et alj 
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; Briiggen 


2003; Dalla Vecchia et alj 120041 


2007 


; Briiggen i 


i Scannapieco 2009|). 



Various authors have implemented self-consistent models of 
black hole growth and AGN feedback in cosmological simulations 
of galaxy groups and clusters (in addition to cooling , star forma- 
tion, an d thermal and chemical feedback from SNe). Spring eTet al.1 
(2005 a ) developed a m odel for quasar mode AGN feedback (see 
also lDi Matteo et al J2005h . w hich was used in cosmolog ical simu- 
lations of 10 galaxy groups bv lBhattacharva et al.l d2008h ■ A model 
for radio mode AGN feedback base d on bubble injection was pro- 
posed by Siiac ki~& Springej j2006t) , which was subsequently ex- 
tended bv lSiiacki et ah (2007) to include quasa r mod e AGN feed- 
back as well. Both Isiiacki & Springell J2006h and Isiiacki et all 
( 2007) performed cosmological simulations of a few massive clus- 
ters with their respective models. 

These studies demonstrated, in a qualitative manner, that AGN 
feedback is effective in reducing the amount of cold baryons and 
star formation in the central regions of groups and clusters. Further- 
more, the gas density is redu ced and the tempera ture is increased, 
elevating the central entropy. Isiiacki etaljj2007l) also showed that 
AGN outflows drive metals from dense, star-forming regions to 
large radii, flattening ICM abundance profiles relative to those pre- 
dicted by a run without AGN feedback. Such trends are precisely 
what is required to reconcile simulations of galaxy clusters with 
observations. 

A more quantitative assessmen t of the impact o f AGN feed- 
back on the ICM was conducted bv lPuchwein et al] J2008h - They 
re simulated a sample of 21 groups and clusters with the scheme 
of lSiiacki et aT] |2007), finding that the model could reproduce the 
observed X-ray luminosity-temperature scaling relation, at least on 
average. However, since their sample size is quite small, it is un- 
clear whether the model can generate a realistic population of CC 
and NCC systems and thus explain the observed scatter about the 
mean relation. In addition, the stellar fraction within the virial radii 

of their simulated objects appears larger than observed. 

Another detailed study was undertaken by iFabian et al.l 

(2010), who resimulated a sample of groups and clusters in a 
cosmologica l settin g, using a model closely related to that of 
Siiack i et alj J2007t) . but with a different implementation of radio 
mode AGN feedback. On group scales, they found that AGN heat- 
ing was able to successfully balance radiative cooling, reproduc- 
ing observed stellar fractions, but the central entropy (at r2soo) 
was about a factor of 2 too high. In addition, their predicted 
group Fe abundance profiles are flat for r > 0.3rsoo, whereas 
observed profiles have a negative gradient out to the largest radii 
for w hich measurements are possible (e.g. iRasmussen & Ponmanl 
2009). There is also an indication that the Fe distribution may be 
too sharply peaked in central regions compared to observations. 
Th e effect of AGN feedback on galaxy groups was also investigated 
by iMcCarthv et al.1 J2010t). who i mplemented the AGN feedback 
scheme of Booth & Schavd l2009h in a cosmological simulation. 
With this model they were able to explain the observed entropy, 
temperature and Fe abundance profiles of groups, as well as ob- 
served X-ray scaling relations. 

For massive clusters, Fabian et al. I d2010h showed that their 
model can reproduce the entropy structure of the ICM, but a fac- 
tor of 3-4 too many stars were formed. The cluster Fe abundance 
profiles they obtained have a shape consistent with that of observed 
profiles, although with a higher n ormalisation, but th e central Fe 
abundance may be over-estimated. [Dubois et all d201 ll) also exam- 
ined the role of AGN feedback in establishing the properties of the 
ICM, using a cosmological AMR simulation of a massive cluster 



with a prescription for jet heating by AGN. The entropy profile of 
their cluster agrees well with that of observed CC clusters if metal- 
cooling is neglected, and when metals are allowed to contribute to 
the radiative cooling, the resulting profile resembles that of a NCC 
cluster instead. However, the metallicity profile of their cluster ap- 
pears steeper than observed. 

1.2 This work 

In this work, we pursue a different, but complementary, approach 
to the theoretical study of galaxy clusters. Instead of undertak- 
ing self-consiste nt hydrodynamica l simu lations, we adopt the hy- 
brid approach of Short & Thoma| d2009l . hereafter SHT09) which 
couples a semi-analytic model (SA model) of galaxy formation to 
a cosmological iV-body/smoothed-particle hydrodynamics (SPH) 
simulation. In this model, the energy imparted to the ICM by SNe 
and AGN is computed from a SA model and injected into the bary- 
onic component of a non-radiative hydrodynamical simulation; see 
SHT09 for details. The main advantage of this approach is that 
feedback is guaranteed to originate from a realistic population of 
galaxies, since SA models are tuned to reproduce the properties of 
observed galaxies. As a conseque nce, the stellar fract ion in massive 
clusters agrees with observations jYoung etalj|201ll) . which is not 
the case in self-consistent hydrodynamical simulations. 

We have extended the mode l of SHT09 to follo w the metal 
enrich ment of the ICM. Note that lCoral d2006l see also lCora et al.1 
2008) have already used a similar hybrid technique to study the 
pollution of intracluster gas by heavy elements. However, they did 
not include energy injection from SNe and AGN, which are likely 
to affect the distribution of metals in the ICM. 

In the model of SHT09, the energy liberated by SN explosions 
and black hole accretion is assumed to be distributed uniformly 
throughout the diffuse gas of the host halo. With this rather ad hoc 
heating model they were able to reproduce observed X-ray scaling 
relations for NCC clusters, but ICM entropy profiles were found 
to be flatter than observed within 0.5 times rsoo (STY10). These 
simulations do not well resolve the core (r < 0.1 rsoo), nor do 
they include radiative cooling that is likely to be important in this 
region, at least for CC clusters. However, we would expect that they 
should be able to provide a much better fit to X-ray observations of 
NCC clusters outside the core. 

The primary goal of this paper is, therefore, to formulate a new 
feedback model that has a clear physical motivation and that is bet- 
ter able to explain the radial variaton of both the thermal and chem- 
ical properties of intracluster gas outside the core of the cluster. To 
help us do this we test a wide variety of different models for SN and 
AGN feedback and metal enrichment, using a selection of X-ray 
data (namely, entropy and metallicity profiles and the luminosity- 
temperature scaling relation) to identify the features that a model 
should possess in order to reproduce the data. 

Our conclusion is that a stochastic heating model, motivated 
by observations of anisotropic AGN outflows, provides a better fit 
to the observed properties of the ICM than more commonplace 
models, such as heating a fixed number of neighbours or heating 
particles by a fixed temperature. Using entirely plausible duty cy- 
cles and opening angles for the jets, it is possible to provide an 
acceptable fit to all available observations with our model. 

Note that the use of SA models means that the feedback is not 
directly coupled to the cooling of the gas - that is why our previous 
work and the bulk of this paper uses non-radiative simulations and 
restricts its attention to NCC clusters. However, towards the end 
of the paper we introduce radiative cooling in an attempt to repro- 
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duced CC clusters. We estimate the degree to which the SA model 
fails to supply the required feedback energy and show that there 
can be a substantial short-fall at high redshift, but that it averages 
to under 10 per cent over the lifetime of the cluster. We are able to 
qualitatively reproduce some CC profiles, but we do not provide a 
detailed quantitative analysis here. 

In this work, we neglect many physical effects such as mag- 
netic fields, cosmic rays, thermal conduction, turbulent mixing, 
etc.. Our principal reason for doing this is to keep the model simple 
and ease interpretation of our results. Some of these may be im- 
portant in the central regions of CC clusters (r < rsoo) but there 
is little evidence that they play a significant role at the larger radii 
that we use to constrain our models. We discuss this further at the 
end of the paper. 

The layout of this paper is as follows. In Section[2] we present 
the details of our hybrid numerical model and describe our cluster 
simulations. We investigate the effect of SN feedback on the ther- 
mal and chemical properties of the ICM in Section [3] and assess 
how our results are affected by different choices of SN feedback 
and metal enrichment models. We show, in agreement with previ- 
ous work, that SN have little impact on the entropy structure of the 
intracluster gas. In Section [4] we examine the impact of additional 
heating from AGN: these can reproduce the correct scaling rela- 
tions but give entropy profiles that are too flat. Our results motivate 
a new, stochastic feedback model based on jet heating, which is 
described in Section [5] In this section, we also discuss what this 
model predicts for the thermal and chemical properties of the ICM, 
and we conduct an exhaustive comparison with observational data 
in Section [6] In Section [7] we demonstrate that our model is capa- 
ble of producing both CC and NCC clusters with the inclusion of 
metal-dependent radiative cooling. Our conclusions are presented 
in Section[8] 

For those readers who are mostly interested in the final model 
itself, rather than the steps used to motivate it, we recommend skip- 
ping Sections[3]and[4] at least on first reading. 



Table 1. The masses, M (in units of h~ 1 Mq), and dynamical tempera- 
tures, fceTdyn (in units of keV), of the 25 clusters used in this study within 
''SOO (second and third columns, respectively), and r2oo (third and fourth 
columns, respectively). Cluster CI is our fiducial cluster, used for most of 
the plots in this paper. 



Cluster name 


-^500 


Tdyn,500 M2OO 


Tdyn,200 


CI 


2.7 x 10 14 


3.9 4.6 x 10 14 


3.9 



C2 


7 


1 X 


10 14 


7 


5 


1 


1 


X 


10 15 


7.1 


C3 


4 


2 x 


10 14 


5 


8 


5 


8 


X 


10 14 


5.4 


C4 


3 


5 X 


10 14 


5 


2 


4 


9 


X 


10 14 


4.8 


C5 


3 


9 x 


10 14 


5 


6 


6 


1 


X 


10 14 


5.3 


C6 


7 


3 x 


10 14 


9 


g 


I 


j 




10 15 


8.6 


C7 


5 


7 x 


10 14 


7 





9 





X 


10 14 


6.5 


C8 


5 


X 


10 14 


5 


7 


7 


1 


X 


10 14 


5.3 


C9 


3 


7 x 


10 14 


5 





5 





X 


10 14 


4.7 


CIO 


3 


9 x 


10 14 


5 


2 


5 


2 


X 


10 14 


4.8 


Cll 


2 


8 x 


10 14 


4 


3 


4 


1 


X 


10 14 


4.0 


C12 


3 


4 x 


10 14 


4 


7 


5 


2 


X 


10 14 


4.5 


C13 


3 


5 x 


10 14 


4 


7 


5 


1 


X 


10 14 


4.5 


C14 


3 


6 x 


10 14 


4 


9 


5 


9 


X 


10 14 


4.6 


C15 


2 


6 x 


10 14 


4 


3 


4 


2 


X 


10 14 


3.9 


C16 


3 


3 x 


10 14 


4 


6 


4 


5 


X 


10 14 


4.3 


C17 


2 


4 x 


10 14 


4 


1 


4 


1 


X 


10 14 


1.0 


C18 


2 


3 x 


10 14 


3 


9 


3 


7 


X 


10 14 


3.5 


C19 


2 


2 x 


10 14 


3 


5 


3 


2 


X 


10 14 


3.2 


C20 


1 


7 x 


10 14 


3 


6 


3 


5 


X 


10 14 


3.4 


C21 


1 


7 x 


10 14 


3 


1 


2 


3 


X 


10 14 


2.8 


C22 


1 


6 x 


10 14 


2 


9 


2 


2 


X 


10 14 


2.6 


C23 


9 


8 x 


10 13 


2 


5 


1 


9 


X 


10 14 


2.3 


C24 


1 


1 X 


10 14 


2 


2 


1 


6 


X 


10 14 


2.1 


C25 


8 


7 x 


10 13 


1 


9 


1 


3 


X 


10 14 


1.8 



We define dynamical temperature as 7d yn = fimn{v 2 ) /3/cb, where 
/imjj ~ 10 — 27 kg is the mean particle mass and (v 2 ) is the mean square 
velocity. 



2 SIMULATIONS 

We make use of hydrodynamical resimulations of a sample of mas- 
sive galaxy clust ers extracted from th e dark-matter-only Millen- 
nium Simulation ( Spr ingel et alj2005bh . Our sample consists of 25 
objects with 9 x 10 13 /i -1 M < M500 < 7 x 10 14 /i _1 M Q and 
forms a subset of the larger sample of 337 groups and clusters res- 
imulated by STY10 for their so-called FO simulation, one of the 
Millennium Gas Simulation^. See STY10 for details of the cluster 
selection procedure. Basic properties of our clusters are listed in 
Tabled 

Following STY10, the feedback model we adopt in our simu- 
lations is the hybrid scheme of SHT09, where a S A model of galaxy 
formation is used to compute the number of stars formed and the 



3 The Millennium Gas Simulations are a series of hydrodynamical simu- 
lations designed to add gas to the dark matter structures found in the Mil- 
lennium Simulation ISpringel et alj|2005bl) . At present, there are 3 simu- 
lations, each of which employs a different physical mechanism for raising 
the entropy of intracluster gas. The first of these is a reference model that 
includes gravitational heating only (the GO run). The second includes ra- 
diative cooling and uniform preheating at z = 4 as a simple model for 
heating from astrophysical sources (the PC ran). The third simulation is the 
FO run, where feedback from galaxies is computed from a SA model using 
the hybrid model of SHT09. 



energy transferred to the ICM by SNe and AGN. We refer the reader 
to STY10 for a full description of the modelling process and simu- 
lation parameters. 

Briefly, we first perform dark-matter-only simulations of each 
region containing a cluster in our sample usin g the massively paral- 
lel TreePM Af-body/SPH code GADGET-2 JSpringell l2005). Viri- 
alised dark matter haloes are identified at each simulation output 
using the friends-of-friends (FOF) algorithm, with a st andard link- 
ing le ngth of 20% of the mean inter-particle separation jDavis et"al] 
1985). Only groups with at least 20 particles are kept, yielding a 
minimum halo mass of 1.7 x lO 10 ^ 1 Mq. Gravitationally bound 
substructures orbiting within these FOF halo es are then found wit h 
a parallel version of the SUBFIND algorithm dSpringel et al.l200ll) . 
From the stored subhalo catalogues we construct dark matter halo 
merger trees by exploiting the fact that each halo will have a unique 
descendant in a hiera rchical scenario of structure formation; see 
Sprin gel et al.l d2005bl) for further details. 

The second stage is to generate galaxy catalogues for 
each resi mulated region by applying the Munich L-Galaxies SA 
model of |Pe Lucia & Blaizo 3 d2007l) to the halo merger trees. 
A full description of the physical proces ses incorporated i n L- 
Galaxies and model param eters is given in lCroton et alj d2006l) and 
|Pe Lucia & Blaizol d2007t) . For each galaxy in these catalogues, we 
use its merger tree to compute the change in stellar mass, A M, , and 
mass accreted by the central black hole, AA/bh, between succes- 
sive model outputs. Knowledge of AM, enables us to incorporate 
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star formation in our simulations as described below. From AM* 
and AA/bh we can also calculate the energy imparted to intraclus- 
ter gas by Type II SNe, AiJsN, and AGN, AEagn, respectively. 
Details are given in SHT09. 

For the purpose of this work, we have extended the model of 
SHT09 to also follow the enrichment of the ICM by metals ejected 
from galaxies in winds. In the L-Galaxies SA model 43% of the 
mass of newly formed stars is instantaneousl y returned, and de - 
posited in the cold gas disc of the host galaxy dCroton et al.ll2006h . 
In other words, the model assumes that metal ejection is instanta- 
neous and does not distinguish between emission from Type II and 
prompt Type la SNe, and that from delayed Type la SNe and AGB 
stars. This will be added in future work. 

In each model galaxy, metals can reside in several distinct 
phases: stars, cold disc gas, hot halo gas and gas ejected by winds 
from the halo into an external 'reservoir'. Only the latter two are 
relevant for the ICM. We define the total mass in metals in diffuse 
gas to be 



M z ,icm = Mz.hot + M z ,, 



(1) 



where Mz.hot and Mz, e j are the mass in metals in hot and ejected 
gas, respectively. 

Once a galaxy falls into a FOF group, becoming a satellite of 
the central galaxy of the halo, all of its metals in hot and ejected gas 
are assumed to be associated with the central galaxy. It follows that 
Afz,icM is non-zero only for central galaxies. Given a particular 
halo at some output redshift z n , we compute the change in metal 
content of the ICM, AMz,icm, since the previous output, z n +i, 
by taking A/z.icm for the central galaxy and subtracting the sum 
of A/z.icm for every galaxy that is a progenitor of any galaxy con- 
tained in the host FOF group and also a central galaxy of a halo at 



Zn + l- 



AM z ,icm = Afz,icM(zn) - ^2 m z,icm(z„+i) 



(2) 



prog. 



The quantity AA/z,icm is used to implement metal enrichment of 
the ICM in our simulations, as described in subsequent sections. 

Finally, we couple the L-Galaxies SA model to hydrodynami- 
cal simulations of our clusters to track the effect of feedback from 
galaxies on the thermal and chemical properties of the ICM. The 
initial conditions for these resimulations are the same as for the 
dark-matter-only runs described above, except that we add gas par- 
ticles with zero gravitational mass. This ensures that the dark mat- 
ter distribution remains undisturbed by the inclusion of baryons, so 
that the halo merger trees used to generate the semi-analytic galaxy 
catalogues will be the same. Gas particles are added at a lower res- 
olution than the dark matter, simply to ease the computational cost 
of our simulations. The resolution we have adopted is sufficient to 
obtain numerically converged estimates of bulk cluster properties 
for systems with T > 2 keV (SHT09). 

Every time an output redshift is reached in our hydrodynam- 
ical simulations, temporary 'galaxy' particles are introduced at 
positions specified by the SA model galaxy catalogue. For each 
galaxy, we know the increase in stellar mass since the last out- 
put, and we remove this mass from the hot phase by converting 
the AiV s tar = AM,/m gas nearest gas particles into collisionless 
star particles, using a stochastic method to ensure that AiV s t a r is an 
integer. Once star formation is complete, we then distribute metals 
and the heat energy available from SNe and AGN amongst neigh- 
bouring gas particles in some way, as described in the following 
sections. Following the injection of metals and entropy, the galaxy 
particles are removed and the simulation continues until the next 



output time, when the process is repeated. The main purpose of 
this paper is to investigate different ways of heating and enrich- 
ing intracluster gas, using X-ray observations of galaxy clusters to 
constrain our models. 

Cluster catalogues are generated at z = from our simula - 
tions using a procedure similar to that of iMuanwong et ail d2002h . 
Full details of our cluster extraction method are given in STY10. 

Following SHT09, we choose to neglect gas cooling processes 
in our hydrodynamical simulations throughout most of this work. 
Although cooling is relatively unimportant for the majority of the 
ICM, we cannot expect to reproduce the low central entropy and 
steep entropy profiles of observed CC clusters, as demonstrated by 
STY 10. However, in Section|7]of this paper, we make a first attempt 
to overcome this limitation of the model of SHT09 by including 
metal-dependent radiative cooling in our simulations. With the ad- 
dition of cooling, we show that it is indeed possible to produce both 
CC and NCC systems using our hybrid approach. 



3 FEEDBACK FROM TYPE II SUPERNOVAE 

In this section we investigate how galactic winds driven by Type 
II SNe shape the chemical and thermal properties of intracluster 
gas. We pay particular attention to how our results are affected by 
varying the feedback scheme in our simulations. Sections |3.1rt3.3| 
describe the models; then in Sections |3 .41 and |331 respectively, we 
show that SNe simply do not provide enough energy to significantly 
alter the entropy and metallicity profiles of the ICM. Most of the 
metals in the ICM originate outside the central cluster galaxy and 
we argue that the metallicity profile in these models is imposed by 
the accretion history of the ICM. 



3.1 Supernova feedback models 

There are two broad classes of SN feedback models deployed 
in numerical simulations: thermal, where the available energy 
is us ed to raise the temperature of neighbouring gas parti- 



cles dKatz 1 19921: iMori et al 
iKav et al.l (20021 : iBrooketai 



1997 



2004 



Thack er & Couchmanl I2OO0I: 
Stinson et alj l200fj) . and ki- 



netic, where neighbouring particles are given a velocity 'kick' 
dNavarro & Whitelll993l; iMihos & Hemquistll 19941: lKawatdl200ll; 
i Kavetal .120 02; Spring e! & Hernquistl2 003; Oppenh eimer & Pavel 
boodlDubois & Tevssieril2008h . 

It is well known that simple thermal feedback schemes fail in 
simulations with cooling since the injected energy is radiated away 
before it has any hydrodynamical effect. This problem is typically 
evaded by suppressing radiative cooling by hand. However, this is 
not an issue for us since cooling processes are not included in any 
of our simulations until Section [7] We have experimented with a 
variety of both thermal and kinetic models, which we now describe. 



3.1.1 Thermal models 

The thermal feedback models employed in our simulations can be 
grouped into three categories, depending on the method used to 
inject the SN energy, AEsn, into the ICM. 

In our first scheme, we simply heat a fixed number of the gas 
particles closest to each galaxy, where the number of neighbours 
heated is iVhcat — 1> 10 or 100. This is the approach typically 
adopted in fully self-consistent hydrodynamical simulations with 
radiative cooling, star formation and thermal SN feedback. 
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The second method we have investigated is to heat all gas par- 
ticles within some sphere centred on each galaxy, where the radius 
of the sphere is assumed to be some fraction, / ra d, of the halo virial 
radius, r2oo- We have explored / ra( j = 0.1, 0.32 and 1, In this 
model, iVheat is defined as the number of gas particles enclosed by 
the sphere. If no neighbours are found, the radius of the sphere is 
increased until a single gas particle is found. 

Our third approach is to heat neighbouring gas particles by a 
multiple, /temp, of the halo virial temperature, Tboo, defined by 

G pm H M 200 



T200 = 



(3) 



2 k B r200 

The values of /temp that we adopt are 1, 3.2 and 10. The number 
of gas particles that can be heated with the available energy is then 

/jw.h(7 - 1)A£ , SN 



N h . 



(4) 



/temp fcBl200?7lgas 

In GADGET-2 the thermodynamic state of each fluid element 
is defined in terms of the entropic function 



(7 - iK 



Pi 



7-1 



(5) 



where m is the thermal energy per unit mass of a particle and pi 
is its density. Supplying heat energy to a gas particle causes Ai 
to increase. Note that A is related to the X-ray gas entropy K via 
K = fj,mji(fj, c mn) J ~ 1 A, where p c mu ~ 1.90 x 10~ 27 kg is the 
mean molecular mass per free electron. 

In each of our three feedback schemes, we heat particles by 
raising their thermal energy by a fixed amount 

AE SN 



Am; = 

■< v heat "Igas 

implemented in GADGET-2 as an entropy boost of 

(7 - l)Aiii 



AAi 



I7-1 



(6) 



(7) 



[max (/ b p200, pi 
The product of the cosmic baryon fraction, /b, and the virial den- 
sity, P200, gives the mean overdensity of baryons within the virial 
radius. If the required -/Vh cat neighbours are not found within a dis- 
tance r2oo of a galaxy and the search radius has to be increased, 
the density of some of these particles may be less than /bP200- By 
using [max (/bp200, Pi)] 7-1 ' rather than p 7_1 , in the denominator 
of equation Q, we are assuming that the amount of energy used 
to heat such particles is Ait;(pi//bP20o) 7-1 < Au^; the rest of 
the energy is taken to be used up as the gas does work expanding 
adiabatically to a density pi < /bP2oo- 

For the first two schemes mentioned above, we have also 
tested an alternative heating model where gas particles are given 
a fixed entropy, rather than energy, boost: 

(7 - 1)AW Aui 



AA t = 



Efir t [ m ax(/bP2oo,P,) 



I7-1 



(8) 



Denser particles close to a galaxy are then heated to a higher tem- 
perature than more distant, lower density particles. 

3.1.2 Kinetic models 

Kinetic SN feedback is implemented in our simulations by assum- 
ing that gas particles closest to a model galaxy are given a velocity 
kick. The number of particles that receive a kick depends on the 
available energy: 

2A£ S n 



iVk 



(9) 



where the wind speed, « w i n( j, is a free parameter. We have consid- 
ered several different values for the wind speed: « W md = 1 km s~ , 
300 km s~ , 600 km s _1 and 1000 km s _ . For any given galaxy, 
we impose the constraint that the wind speed cannot be less than 
the virial speecQ of the host halo, i>2oo , so that the case u w ind = 1 
km s _1 is equivalent to assuming that material is ejected at the 
virial speed. To ensure that iVkick is an integer, we draw a random 
number r uniformly from the unit interval and compare it with the 
fractional part of iVkick : if r is less (greater) than the fractional part 
of iVkick, we round iVkick up (down) to the nearest integer. 

The velocity of each kicked particle is modified according to 



v + v wind n, 



(10) 



where n is a unit vector that is either oriented in a random direction 
on the unit sphere, or in the direction from the galaxy to the wind 
particle. 

We have only studied kinetic feedback models where the 
wind speed is a constant for all galaxies. Similar models 
are often employed in self-consis t ent hydrodynamical simula- 
tions (e.g. lNavarro& White! 1 19931; ISpringel & Hernquistl 120031: 
iDalla Vecchia & Schavd 12008 ). However, there are other possi- 
bilities, such as momentum-driven winds, where t he wind speed 
scales with the galaxy ve locity dispersion (e.g. Ilvlartinl 120051 : 
lOppenheimer & Dave1l2006t) . and models where the outflow veloc- 
ity increases with galactocentric radius dSteidel et alfcOld) . 



3.2 Metal enrichment models 

We distribute metals amongst gas particles in our simulations as 
follows. For each model galaxy, all gas particles contained within a 
sphere centred on the galaxy are identified. As before, the radius of 
the sphere is chosen to be a fraction, /z,rad, of the halo virial radius, 
where /z,rad = 0.1, 0.32 or 1. The metals in diffuse gas produced 
by the galaxy are then shared evenly amongst these particles, so 
that the metal mass associated with each particle, mz,i, increases 
by an amount 

AMz.icm 



Amz = 



(11) 



iV snrich 

where iVemich is the number of gas particles inside the sphere. 

Given the total mass in metals for a gas particle, we could then 
define its metallicity simply as 

mz,i 



-'part,! 



(12) 



which we refer to as the particle metallicity. However, in this work , 
we prefer to use the smoothed metallicity JOkamoto etd]|2005t 
iTornatore etal]|2007r) , defined by 



pz,i 
Pi 



(13) 



where the smoothed metal mass density, pz,i, is computed in an 
analogous way to the standard SPH density estimate: 



pz.i = mz,jW(\rj - rj\, hj 
3=1 



(14) 



Here iV S ph = 64 is the number of SPH smoothing neighbours and 
W is a spherically-symmetric smoothing kernel, which depends 

4 We define the virial speed of a halo to be the circular velocity at the virial 
radius. 
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Figure 1. Entropy profiles for cluster CI resimulated with different imple- 
mentations of supernova feedback (coloured lines; see legend for model 
details). Note that the profile obtained from a gravitational heating only 
model (GO) is also shown. For comparison, we display observed profiles of 
similar-mass CC (dashed grey lines) and NCC (solid black lines) clusters in 
the REXCESS sample (PAP10). 



Figure 2. Emission- weighted metallicity profiles for cluster C I resimulated 
with the same kinetic supernova feedback model, but varying the radius 
of the region within which metals are injected (solid coloured lines). See 
the legend for details of the metal enrichment models adopted. Observed 
profiles of CC (dashed grey lines) and NCC (solid black lines) clusters from 
the sample of MATH are also shown. 



upon the separation of particles i and j, |r< — r 3 ■ |, and the smooth- 
ing length of particle i, hi. The smoothed metallicity of a gas par- 
ticle is updated whenever its SPH density is calculated. Once a gas 
particle is converted into a star particle, its smoothed metallicity re- 
mains fixed for the rest of the simulation. Note that the metallicity 
of particles does not affect the gas dynamics in our non-radiative 
simulations. 

3.3 Naming conventions 

Table [2] lists all 23 of our SN feedback models. Note that each 
model, including the reference gravitational heating only model 
(GO), follows the conversion of gas into stars as dictated by the 
underlying SA model. 

3.4 Entropy profiles 

To test the effect of different implementations of SN feedback on 
the entropy structure of the ICM, we have resimulated our fiducial 
cluster, CI, with each of our models. Figure [T] shows the result- 
ing entropy profiles. For comparison, we also show the observed 
entropy profiles o f CC and NCC clusters in the REXCESS sample 
dPratt et aljioicl hereafter PAP 10) . To facilitate a fair comparison, 
we only plot profiles of observed clusters that have a mass, M500, 
within 20% of that of cluster CI. 

The main point to note is that all of our models yield almost 
identical entropy profiles that are in good agreement with the pro- 
file obtained from the reference GO run. In all cases, the profiles 
scale approximately as K oc r 1 2 for r > O.lrsoo, consi stent with 
spherical accretion models (e.g. iTozzi & Normartl200lh and cos- 
mological simu l ations that includ e gravitational heating only (e.g. 
IVoit et al.ll200llNagai et alj|2007l) . For r < O.lrsoo, the entropy 
profiles flatten off significantly, exhibiting a small spread in central 
entropy. 

Compared to the observed entropy profiles of CC clusters, 



the profiles predicted by our models have a steeper slope at r > 
O.lrsoo, and the normalisation is systematically too low. In the case 
of NCC clusters, it is evident that none of our models can explain 
the shallow profiles characteristic of these systems. 

We have checked that these results hold for other clusters in 
our sample, so we conclude that SNe have a negligible impact 
on the thermodynamical properties of intracluster gas and, further- 
more, the manner in which the feedback energy is injected is unim- 
portant. 



3.5 Metallicity profiles 

The metal enrichment model has only one free parameter, /z.rad, 
which controls the radius of the spherical region about a galaxy 
in which metals are injected. We have resimulated cluster CI with 
three different values of /z, ra d, fixing the SN feedback scheme to 
be the kinetic model where gas particles are kicked in a random 
direction with « W md = 600 km s _1 . Figure|2]shows the emission- 
weighted metallicity profiles that result, along with observed Fe 
abundance profiles of CC and NCC clusters from Mats ushital d20 1 1[ 
hereafter MATH). We plot all observed clusters with a mass above 
80% of that of cluster CI, in order to obtain a reasonable number 
of both CC and NCC objects (most NCC clusters in the sample of 
MATH are considerably more massive than cluster CI). 

We note that it is difficult to directly compare the metallic- 
ity profiles of our clusters with those of observed clusters from 
MATH, for several reasons. Firstly, the observed profiles are Fe 
abundance profiles, but our simple metal enrichment model does 
not include the contribution from Type la SNe, a major source 
of Fe, nor does it track the production of individual chemical el- 
ements. Secondly, we could, in principle, adjust the yield in the SA 
model underpinning our simulations, which would allow us to al- 
ter the normalisation of our metallicity profiles. For these reasons, 
we focus on the shape of metallicity profiles, instead of their nor- 
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Table 2. Supernova feedback models. Unless otherwise stated, the radius for metal injection is r200 (/z.rad = 
1). Note that all models, including the GO model, follow the conversion of gas into stars. 



Model name 


Type 


Energy injection method 


Comments 


GO 






Gravitational heating only 


SN_Th_N7V hcat 

SN_Th_R/ rad 

SN.Th.T/ tcmp 


Thermal 
Thermal 
Thermal 


Fixed energy 
Fixed energy 
Fixed energy 


AW = 1, 10, 100 
/ rad = 0.1,0.32, 1 
/temp = 1, 3.2, 10 


SN_En_NAW 
SN_En_R/ rad 


Thermal 
Thermal 


Fixed entropy 
Fixed entropy 


AW = 1, 10, 100 
/ rad = 0.1,0.32, 1 


SN_KiR_Vt> wind 
SN_KiD_Vu w ind 


Kinetic 
Kinetic 


Velocity kick, random 
Velocity kick, directed 


Dwind/kms" 1 = 1, 300, 600, 1000 
iiwind/kms -1 = 600 


SN_KiR_Z/ z , rad 


Kinetic 


Velocity kick, random 


/z,rad = 0.1, 0.32, 1; ^jnd/kms- 1 = 600 



malisation, when assessing the impact of different feedback and 
enrichment models on the ICM enrichment pattern. 

It is apparent from Figure [2] that varying the metal injection 
radius has a large impact on ICM metallicity profiles in core re- 
gions, r < 0.2ri8o. If metals are injected in a concentrated fashion 
(/z.rad = 0.1), then we see a sharp peak in the metal distribution 
within that region that is not reflected in the observational data. As 
the injection radius is increased, the gradient of the profile becomes 
progressively shallower until, when /z.rad = 1, the slope provides 
a good match to that of the profiles of observed NCC clusters, ex- 
cept the few systems that have a high central abundance more typi- 
cally found in CC clusters. It is possible that these systems are CC 
remnants. Recall that radiative cooling is not included in our simu- 
lations so we do not expect to be able to reproduce the abundance 
peaks seen in the core regions of CC clusters. 

A metal enrichment model where metals are distributed 
throughout the halo can be justified if the bulk of the metals found 
in intracluster gas were brought in by infalling material, rather than 
being produced by star formation in the central galaxy of the halo. 
To check whether this is the case, we have modified the L-Galaxies 
SA model to follow what fraction of metals in diffuse halo gas are 
produced by the central galaxy of the halo. Figure[3]shows this frac- 
tion as a function of halo virial mass for all 25 clusters in our sam- 
ple, and for a selection of halos with M200 ^ 10 n /i -1 M© taken 
from the Millennium Simulation galaxy catalogues. For all of our 
clusters, the fraction of metals in hot halo gas produced by the cen- 
tral galaxy is less than 5%, implying that nearly all of the metals in 
diffuse gas are indeed accreted. Note that this may change some- 
what when we extend L-Galaxies to track the time-dependence of 
metals returned by Type la SNe and AGB stars, as some of the 
metal production will be delayed until after the formation of the 
central cluster galaxy. 

Ideally one would like to inject metals locally about satel- 
lite galaxies falling into a halo, rather than distributing them uni- 
formly throughout the halo . However, this is not possible with the 
|Pe Lucia & Blaizotl J2007h version of L-Galaxies since all the met- 
als in the diffuse gas associated with a galaxy are assumed to be in- 
stantaneously stripped once it becomes a satellite galaxy. In future 
work we plan to switch to a different treatment of satellite galaxies 
whereby hot gas is gradually remove d from infalling galaxies b y 
tidal and ram-pressure stripping (e.g. iHenriq ues & Thomas feoiOl) . 

We now examine whether changing the SN feedback scheme 
affects ICM metallicity profiles. To do this, we have resimulated 
our fiducial cluster with all of the models described in Section[3Tl 
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Figure 3. Fraction of metals in hot halo gas produced by the central galaxy 
of the halo, /z.ccnt. as a function of virial mass. The black points are for 
a subsample of halos with M200 ^ 10 11 /! - 1 Mq extracted from the Mil- 
lennium Simulation galaxy catalogues. The red circles correspond to our 25 
resimulated clusters. For these massive systems, under 5% of metals in the 
hot gas are produced by the central galaxy, implying that nearly all of the 
metals are accreted. 



fixing /z,rad = 1 in each case. Figure [4] compares the emission- 
weighted metallicity profiles obtained from our cluster simulations 
with those of the same observed clusters from the MATH sample. 

The main point to note is that the metallicity profiles obtained 
from all our various simulations are essentially the same and we 
have checked that this conclusion remains valid when metals are 
injected in a concentrated manner (/z,rad = 0.1). It follows that 
SN feedback has no impact on the metal distribution in clusters, 
and the precise way in which the energy available from SNe is used 
to heat intracluster gas is irrelevant. 

The metallicity profiles of clus ters in the presence of super- 
nova feedback were investigated by iTornatore et al.l |2007). They 
also found that changing the SNe feedback rate makes little differ- 
ence to the slope of the profiles (although it does change the nor- 
malisation). We show in FigureQT]and Section |531 that the stronger 
AGN jet feedback can have a larger effect. 
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Figure 4. Emission-weighted metallicity profiles for cluster CI resimu- 
lated with different supernova feedback schemes, assuming that metals are 
distributed uniformly throughout the halo (coloured lines; see the legend 
for feedback model details). We also show profiles of CC (dashed grey 
lines) and NCC (solid black lines) clusters from the observational sample 
of MATH. 

3.6 Summary 

Our study so far has revealed that feedback from SNe has a negligi- 
ble effect on both ICM entropy and metallicity profiles, regardless 
of the manner in which the energy is assumed to be transferred to 
the gas. 

In light of this freedom, we choose our fiducial SN feedback 
scheme to be the kinetic model SN_KiR_V600, where gas particles 
are given a kick in a random d irection with i> W md = 600 km s ~ . 
This is similar to the model of iDalla Vecchia & Schavd j2008l) . A 
wind speed of 600 km s~ is consistent with o bservations of loca l 
(e.g. IVeilleux et ai]|2005t) and z ~ 2-3 (e.g. ISteidel et alj|201fj|) 
starburst galaxies. 

The metallicity profiles reflect the manner in which metals 
are injected into the diffuse gas. For our fiducial metal enrich- 
ment model we assume that the metals ejected from galaxies are 
distributed uniformly throughout the entire halo, since this gives a 
good match to the slope of the metallicity profiles of observed NCC 
clusters. This model is justified by the fact that nearly all of the met- 
als in intracluster gas are accreted, rather than being produced by 
the central galaxy of the halo. 

The model that forms the basis for the rest of the work pre- 
sented in this paper is thus SN_KiR_Zl. 



4 FEEDBACK FROM ACTIVE GALACTIC NUCLEI 

As shown in the previous section, the heating of intracluster gas 
by stellar feedback alone clearly cannot account for the excess en- 
tropy observed in cluster cores, indicating that an additional feed- 
back mechanism must be at play. The favoured candidate is the 
energy liberated by the accretion of gas onto central supermassive 
black holes at the centres of galaxies. Our goal in this section is to 
assess how the properties of the ICM are altered by the inclusion of 
this extra heating from AGN, and how our results are affected by 
different numerical implementations of AGN feedback. 



In Section l4~71 we describe simple AGN heating models, sim- 
ilar to those found in the literature, then in Section 14.21 we use 
comparisons with observed entropy profiles to conclude that none 
of these models are entirely satisfactory. One model with extreme 
wind speeds does provide an adequate fit to the data and that moti- 
vates the new stochastic heating model developed in Section|5] 

4.1 AGN feedback models 

The amount of energy available from AGN heating, AEagn, is 
not arbitrary but is set by the model described in Secton 3.1.2 
of SHT09. Although that paper considered only a single heating 
model, we find that the global X-ray lumnosity-temperature rela- 
tion is dependent mainly upon the normalisation of the heating, 
and is relatively unaffected by the particular manner in which the 
heat is injected. That can have a large effect on the entropy profiles, 
however, as we show below. 

4.1.1 Thermal models 

The first set of thermal AGN feedback models that we have tested 
are identical to the thermal SN feedback models described previ- 
ously in Section [3. 1.11 except with A£sn replaced by A£agn- 
Similar prescriptions for A GN feedback have been employed in 
numerous other works (e.g. [Springel et al]|2(X)5al : IdT Matteo et al.l 
l2008l:lBooth & Schavej2009l : lFabian et aU2010h . 

4.1.2 Kinetic models 

We implement kinetic AGN feedback in our simulations in the 
same way as kinetic SN feedback; see Section 13.1.21 The only 
differences are that the number of particles kicked (equation [9} 
now depends on the energy available from black hole accretion, 
A£agn, rather than that available from SN explosions, AEsn, 
and we have adopted larger wind speed values, u w i n( j = 1000 
km s _1 , 4500 km s _1 a nd 20000 km s -1 , i n line with measured 
AGN outflow velociti e s jPounds et alj|2003l; IChartas et alj|2003l; 
Crenshaw etal] 120031 ; iGangulv & Brothertonl 120081 ; iDunn et al.l 

2o"ioh . 

4.1.3 Naming conventions 

Table|3]lists all of our various AGN feedback models. In each case 
the SN feedback and metal injection schemes are the same as for 
model SNJCiRJZl. 

4.2 Entropy profiles 

In order to assess how sensitive the thermodynamical properties of 
the ICM are to different implementations of AGN feedback, we 
have resimulated our fiducial cluster with each of our 15 thermal 
and 4 kinetic AGN feedback models. 

4.2.1 Thermal models 

The entropy profiles obtained from our thermal models are dis- 
played in Figure [5] For comparison, we also show the profile pre- 
dicted by our fiducial SN feedback model (SN_KiR_Zl), and en- 
tropy profiles of observed clusters of a similar mass. 

As for the SN feedback, all of the thermal models give very 
similar results, even though there are large differences in the way 
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Table 3. AGN feedback models. In each case, supernova feedback is implemented using a kinetic model where 
particles neighbouring a galaxy are given a kick in a random direction with velocity 600 km s — , and ejected 
metals are assumed to be distributed uniformly throughout the entire host halo. 



Model name 


Type 


Energy injection method 


Comments 


AGN.Th.NAW 
AGN.Th.R/ rad 
AGN.Th_T/ t cmp 


Thermal 
Thermal 
Thermal 


Fixed energy 
Fixed energy 
Fixed energy 


AW = 1, 10, 100 
/ rad = 0.1,0.32, 1 
/temp = 1, 3.2, 10 


AGN.En.NAW 
AGN.En.R/ rad 


Thermal 
Thermal 


Fixed entropy 
Fixed entropy 


AW = 1, 10, 100 
/ rad = 0.1,0.32, 1 


AGN.KiR.ViWd 
AGN.KiD_VtWind 


Kinetic 
Kinetic 


Velocity kick, random 
Velocity kick, directed 


^wind/fans" 1 = 1000, 4500, 20 000 
Dwind/kms- 1 =20000 




Figure 5. Entropy profiles for cluster CI resimulated with different imple- 
mentations of thermal AGN feedback (coloured lines; see the legend for 
model details). For comparsion, we also show the profile obtained from a 
ran with kinetic supernova feedback only, model SN_KiR_Zl, and observed 
profiles of CC (dashed grey lines) and NCC (solid black lines) clusters in 
the REXCESS sample (PAP10). 



Figure 6. Entropy profiles for cluster CI resimulated with different kinetic 
AGN feedback models (coloured lines; see legend for model details). For 
comparsion, we also show the profile obtained from a ran with kinetic su- 
pernova feedback only, model SN_KiR_Zl, and observed profiles of CC 
(dashed grey lines) and NCC (solid black lines) clusters in the REXCESS 
sample (PAP10). 



in which the available energy is shared amongst gas particles in 
the various schemes. Essentially, the profiles follow the predicted 
r i. 1—1.2 sca j m g at i ar g e ra (Jii, but as we move in towards the core 
they begin to flatten off at r ~ 0.5 — 0.6rsoo- At radii interior to 
this, the slope is shallower than seen in either CC or NCC clusters, 
leading to an over-estimate of the central entropy. Simple preheat- 
ing models predict similarly large isentropic cores at z = (e.g. 
STY 10). 

4.2.2 Kinetic models 

Figure [6] compares the entropy profiles predicted by our kinetic 
AGN feedback models with the same set of observed cluster pro- 
files as in Figure[5] 

In cluster outskirts, r ~ rsoo, the models give similar results, 
but there are clear difference at radii less than this. For the low- 
est wind speed, t> W md = 1000 km s -1 , we see a very flat entropy 
profile, with a hint of an entropy inversion in the core. As the wind 
speed is increased, the entropy profile steadily steepens, providing a 
good match to observed NCC cluster profiles when n W md = 20 000 



km s" 1 . Almost identical results are obtained when kicks are im- 
posed in the direction from the galaxy to the wind particle, rather 
than in a random direction. 

To understand this behaviour, we have checked how the tra- 
jectories of kicked particles are affected by variations in the wind 
speed. This is done by identifying the main progenitor of our cluster 
at high-redshift using the halo merger trees, selecting all gas parti- 
cles within rsoo of this object that have just received a kick, then 
tracking the cluster-centric positions of these particles to z = 0. 
For a high wind speed, f w i nd = 20 000 km s _1 , the available AGN 
energy is only sufficient to kick a small number of particles and we 
find that their large momentum carries them beyond rsoo- As we 
reduce the wind speed to u w i nd = 1000 km s" 1 , the number of 
particles kicked increases but their momentum gain is smaller, so 
they do not escape from the cluster core before their kinetic energy 
is converted to thermal energy. This leads to an increase of the gas 
entropy in the central regions, establishing a flat entropy profile as 
seen in Figure [6] 
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4.3 Summary 

Several interesting results have emerged from our study of the ef- 
fect of different AGN feedback models on the thermal properties of 
the ICM. 

We have found that simple thermal feedback schemes, based 
on heating a fixed number of particles, heating particles within a 
fixed fraction of the virial radius, or heating particles by a fixed 
fraction of the virial temperature, all heat the gas in cluster central 
regions excessively, leading to a higher core entropy than observed. 
All 15 of the thermal models we have tested give very similar re- 
sults. 

When AGN feedback is implemented in a kinetic manner, 
ICM entropy profiles are found to be sensitive to the wind speed 
adopted. For low wind speeds, the resulting entropy profiles are too 
flat, as in the thermal case. This is because the available energy is 
shared amongst a large number of particles, and kicked particles 
do not have sufficient momentum to escape central cluster regions 
before their kinetic energy is thermalised. 

As the wind speed is increased, the number of particles kicked 
decreases as 1 /v% ind and kicked particles are able to reach larger 
cluster-centric radii before thermalisation of their kinetic energy. 
Consequently, more low-entropy material remains in core regions, 
so entropy profiles become progressively steeper, approaching ob- 
served ones. For u w i n d = 20 000 km s , the predicted profiles 
agree well with observed profiles of NCC clusters. However, such 
a high wind speed is perhaps physically unrealistic. 

From our discussion, it seems that the key ingredient of a suc- 
cessful AGN feedback model must be to ensure that only a small 
fraction of particles in central cluster regions are heated/kicked, 
so that these particles have sufficient entropy/momentum to reach 
cluster outskirts, leaving low-entropy metal-rich gas behind in the 
core. In the next section we formulate a new AGN feedback pre- 
scription that has this desired feature, and is motivated by the ob- 
served interaction of AGN with their environment. 



5 A NEW MODEL FOR FEEDBACK FROM ACTIVE 
GALACTIC NUCLEI 

There is a growing body of observational evidence that AGN 
feedback may be mostly related to radio-loud AGN. In the local 
universe, observations of galaxy groups and clusters often show 
X-ray cavities coincident with l obes of radio emiss i on linked to 
the central galaxy by radio jet s dBlanton et alj|200ll ; 



McNamara et alj 120051: iFabian et alj 120061: 



Birzan etal 



, Morita et al. 

Jetha e t al. 2008; Gastalde lloet alj|2009l ; iDong et alj|2010l 



2004 
2006 

Giacintucci ct al .201 l|). It is thought that these bubbles are inflated 
by the central AGN, and may provide an efficient means of remov- 
ing cool, enriched gas from cluster cores as they rise buoyantly 
through the cluster atmosphere, thus quenching star formation. 

At high redshift, z ~ 2 — 3, emission-line kinematics of ra- 
dio galaxies based on rest-frame optical integral-field spectroscopy 
have revealed powerful bipolar outflows with kinetic energies 
equivalent to 0. 2% of the rest mass of the central supermassive 
black hole (e.g. iNesvadba et alj|200^. l2008h . These AGN-driven 
winds are energetic enough to remove copious amounts of gas 
from the host galaxy, preventing further accretion onto the black 
hole and suppressing star formation. Large-scale energetic outflows 
have also been observe d in z w 2 ultra-luminous infrared galaxies 
dAlexander et al.l2010h . a galaxy population potentially an order of 
magnitude more common than distant radio galaxies. 



Although it is not yet fully understood how the energy re- 
leased by black hole accretion is transferred to the surrounding gas, 
the observational data suggests that the energy is input in a direc- 
tional manner, via jets or collimated outflows, rather than isotrop- 
ically. To reflect this, we have developed an anisotropic, stochas- 
tic heating model where only some of the gas particles neighbour- 
ing a galaxy are heated per duty cycle of the AGN. We note that 
higher-resolution models of feedbac k from AGN in clu ster cores 
also favour anisotropic heating (e.g. iGaspari et al.ll2012h . 

In Section [5TT1 we describe this heating model in detail, then 
in Section [5721 we use observed entropy profiles and X-ray scaling 
relations to determine optimal model parameters. 



5.1 Stochastic AGN feedback model 

The basis of our new model for AGN heating is as follows. For each 
galaxy, we first identify all gas particles contained within a sphere 
centred on the galaxy, where the radius of the sphere is some frac- 
tion, / ra d, of the halo virial radius. We then assume that the proba- 
bility that any of these particles has been heated by AGN feedback 
during the time elapsed, At, since the previous SA model output is 

\At/t d 



i-(i-/dut y r t/tduty , 



(15) 



where /duty is a parameter controlling the fraction of particles 
heated over the AGN duty cycle, tdutv Based on observational 
data, we take tdutv = 1 8 yrs (e.g. lBirzan et alj2 004; Fabia rTetal] 
l2004IJetha et al.12008 "). With our choice of SA model output times 
we then have 2 < At/t^uty < 4 for z < 3. 

For each gas particle neighbour, we draw a random number 
r uniformly from the unit interval and compare it with Phcat: if 
r < Pheat the particle is given an entropy boost 



( 7 - 1)AE A 



Pheat XL=r* [ max (/bP200,Pj 



17-I ' 



(16) 



and if r ^ P nca t the particle is not heated. By including the heating 
probability Ph cat in the denominator of equation <16t , we ensure 
that the total amount of energy injected into the gas is (approxi- 
mately) the same for different choices of /duty 

We have also experimented with supplying the AGN heat en- 
ergy to particles as a fixed energy boost. However, this makes virtu- 
ally no difference to our results so we do not discuss these models 
hereafter. 

There are two free parameters in our model: / ra d and /duty 
The values of these parameters we have tested in this work are 
/ rad = 0.1, 0.32 and 1, and f duty = 10~ 4 , 10~ 3 , 10~ 2 and 
10 _1 . Note that in the case /duty = 1 our model reduces to 
AGN.Th.R/ rad . 

It is interesting to link the parameter /duty to the opening 
angle of AGN jets. If we make the simple approximation that 
large-scale AGN outflows can be treated as biconical jets, each 
with opening angle 28, then it follows that cos# = 1 — /duty 
For the range of values of /duty tested here, this corresponds to 
1° < 6 < 26°. 

Table [3] lists all of our stochastic AGN feedback models. To 
distinguish these from the AGN heating models of the previous sec- 
tion, we have given them the label JET. 



5.2 Entropy profiles and scaling relations 

We now investigate whether our new physically-motivated stochas- 
tic AGN feedback scheme yields a better match to observed cluster 
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Table 4. Stochastic AGN feedback models. In each case, supernova feedback is implemented using a kinetic 
model where particles neighbouring a galaxy are given a kick in a random direction with velocity 600 km 
s _1 . Unless otherwise stated, the radius for energy and metal injection is 7*200 (/rad = 1 ar, d /z rad = 1> 
respectively), and the fraction of particles heated per AGN duty cycle is /duty = 10 — 2 . 



Model name 


Type 


Energy injection method 


Comments 


JET_R/ rad _D/duty 


Stochastic 


Fixed entropy 


/ r ad =0.1,0.32, 1 








/duty = 10- 4 , lO" 3 , 10- 2 , 10- 1 


JETJZ/z.rad 


Stochastic 


Fixed entropy 


/Z,rad =0.1,0.32, 1 



E 

o 
> 




0.01 0.10 1.00 

r / r 500 



Figure 7. Entropy profiles for cluster CI resimulated with our stochas- 
tic AGN feedback model for different values of the parameter / rad (solid 
coloured lines; see legend for model details). The other parameter in the 
model, /duty. i s fixed at 10 -2 . The profiles of observed CC (dashed 
grey lines) and NCC (solid black lines) clusters in the REXCESS sample 
(PAP 10) are also displayed for comparison. 

profiles than the simple thermal and kinetic models discussed in the 
previous section. 

Recall that our stochastic model has two free parameters: / ra d, 
which governs the radius about a galaxy in which energy is in- 
jected, and the fraction, /duty, of neighbouring gas particles that 
are heated per AGN duty cycle. The first issue to address is how 
varying these parameters affects cluster properties. We then iden- 
tify an optimal choice for these parameters by using a selection of 
observational data to constrain our model. 
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Figure 8. The X-ray luminosity-temperature scaling relations predicted by 
our stochastic AGN feedback model with / rad = 0.1 (asterisks), 0.32 (tri- 
angles) and 1 (circles), keeping the other model parameter, /duty fixed 
at 10 — 2 . See the legend for model names. X-ray properties are computed 
within r,5oo- For comparative purposes, we also plot observational data 
for CC (diamonds) and NCC (squares) clusters in the REXCESS sample 
(PCA09). 

excessively large and excessively small values of / ra d will both 
leave behind low-entropy core particles. However, the variation is 
small and each of the chosen values yields an adequate match to 
the observed relation for NCC clusters, with / ra d = 1 providing 
the best match of the three. 



5.2.1 The effect of changing / ra d 

Figure [7] shows the effect of varying / ra d on the entropy profile of 
cluster CI. In each case, / du t y is kept fixed at 1CP 2 . It is apparent 
that the entropy structure of the ICM is relatively insensitive to the 
choice of / ra d, with only small differences between the three dif- 
ferent runs. The best match to observed NCC cluster profiles arises 
when / ra( j = 1, in which case we find excellent agreement with the 
observational data. 

Figure[8]shows the Lx-Tsi relation for our full 25-cluster sam- 
ple for the three different values of / ra d tested. Here also, we can 
see that the three different models predict a very similar Ly^-T B \ re- 
lation. Note that the trend in Lx-Tsi with / ra d is not monotonic: 



5.2.2 The effect of changing /duty 

We now turn our attention to the effect of the parameter /duty We 
have done four runs, with /duty = 10~ 4 , 10~ 3 , 10~ 2 and 10 -1 , re- 
spectively, keeping / ra d fixed at unity. The entropy profile of clus- 
ter CI in each case is displayed in Figure|9] It is immediately clear 
that varying /duty has a much larger effect on the entropy of intra- 
cluster gas than / ra d- As /duty is increased from 10~ 4 to 10" 1 , the 
slope of the entropy profile at radii r < 0.4rsoo becomes progres- 
sively shallower. For /duty = 1CP 4 the slope is too steep compared 
to that of observed NCC cluster profiles, whereas it is too flat for 
/duty = 10" \ The values / duty = 10" 3 and / du t y = 10" 2 both 
give a good match to the observational data for NCC clusters. 
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Figure 9. Entropy profiles for cluster CI resimulated with our stochastic 
AGN feedback model for several choices of the parameter /duty (solid 
coloured lines; see legend for model details). The other model parameter, 
/rad, is set to unity. For comparison, we also show the profiles of observed 
CC (dashed grey lines) and NCC (solid black lines) clusters in the REX- 
CESS sample (PAP10). 

To explain the variation in cluster entropy profiles with /duty, 
we have again examined what happens to particles that are heated 
by AGN feedback in each of our runs. In the case where /duty = 
10 -4 , the probability of a particle being heated is low, but any parti- 
cle that is heated receives a large entropy boost since P nC at appears 
in the denominator of equation U6b . The high entropy of heated 
particles causes them to rise buoyantly to large cluster-centric radii, 
T > ^500, leaving the entropy profile in the core relatively undis- 
turbed compared to a run with SN feedback only. When /duty is in- 
creased to 10" , many more particles in central regions are heated 
by AGN feedback since Phcat is larger, and the entropy boost they 
are given is smaller. Accordingly, the distance they move outwards 
from the core is less, resulting in a higher central entropy and a 
flatter profile. 

The large impact of /duty on cluster entropy profiles is re- 
flected in the L^-T B i scaling relation. This is demonstrated in Fig- 
urellOlwhere we show the Lx-T s i relation for our full cluster sam- 
ple predicted by each of our four models. For /duty = 10 -4 , the 
low central entropy causes an enhanced X-ray luminosity, so all of 
our simulated clusters lie well above the mean observed relation for 
NCC clusters. In fact, the predicted relation in this case resembles 
the observed relation for CC clusters, although this is artificial since 
we have not included cooling processes. As /duty is increased, the 
normalisation of the Lx-T s i relation decreases and the slope be- 
comes steeper. A good match to the observed NCC cluster Lx-T B \ 
relation is obtained when /duty = 10 -2 . For larger values of /duty, 
the relation is too steep, so that low-temperature systems have too 
low a luminosity for their mass. This is because the AGN heating 
has raised the core entropy in these systems to an excessive level. 

5.2.3 Identifying optimal parameter values using observations 

The next issue to address is whether observational data can help us 
to constrain the two free parameters of our stochastic AGN feed- 
back model. For this analysis, we have resimulated all 25 clusters 
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Figure 10. The X-ray luminosity-temperature scaling relations predicted by 
our stochastic AGN feedback model with /duty = 10~ 4 (asterisks), 10~ 3 
(triangles), 10 — 2 (crosses) and 10" 1 (circles), keeping the other model pa- 
rameter, / ra d, set to unity. See the legend for model names. X-ray properties 
are computed within rsoo- Observed CC (diamonds) and NCC (squares) 
clusters from the REXCESS (PCA09) are also displayed. 



in our sample (C1-C25) using different combinations of the pair of 
parameters (/ ra d, /duty)- The values adopted are / ra d = 0.1, 0.32 
and 1, and / duty = 10" 4 , 10" 3 , 10" 2 and 10"\ giving a grid of 
12 models in total. 

We assess the suitability of each model by testing how well 
it reproduces the observed scaling of three fundamental ICM ob- 
servables with spectroscopic temperature: (i) the entropy profile 
normalisation, which we take to be the entropy at riooo (typically 
about 0.7r5oo), (ii) the entropy profile shape, defined as the ratio 
of the entropy at riooo to the entropy measured at 0.1r2oo (i.e. the 
central entropy), and (iii) the X-ray luminosity. Again, the source 
of the observational data is the REXCESS. Note that we are only 
aiming to match the observed scaling relations for NCC clusters 
since cooling is not included in our simulations at this stage. 

In the following analysis, we neglect any clusters in our simu- 
lated sample that have large amounts of substructure. To differenti- 
ate between dynamically relaxed and disturbed systems, we use the 
substructure statistic 

S = l Xc ° m ~ Xc l [ ( i7) 

f500 

where x c is the location of the dark matter potential minimum and 
Xcom is the cen tre of mass of the cluster within rsoo- Following 
lKavetal.l ( l2007T) . we say that a cluster is disturbed if S > 0.1, and 
relaxed otherwise. 

For each scaling relation, the criterion we use to test how well 
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our models reproduce the mean observed relation is the \ 2 statistic 



Table 5. x test probability values for the Lx-T s i scaling relation. 



4- J2 i lo sio \E(zTYr] 



(18) 



/j-isim 

" QlOSl ° I 5teV I ~ l "^" r '" <" ■ 

where N^ lm is the number of (relaxed) simulated clusters and Y = 
K(riooo), K(nooo)/ K(0.1r2oo) or Lx, depending on the relation 
being considered. The quantities Co and a are the normalisation 
and slope of a power-law fit to the corresponding observed relation, 



E{z) n Y° hs = Co 



Tiobs 
X 

5keV 



(19) 



obtained by using the BCE S orthogonal linear regression method 
dAkritas & Bershadvlll996l) in log-log space, taking into account 
the errors in both T x bs any y obs . The normalisation Co has units of 
h~ 1/:i keV cm 2 and 10 44 /i -2 erg s" 1 for Y = if(riooo) and L x , 
respectively, and is dimensionless for Y = -R"(riooo)/-ff(0.1r2oo)- 
The factor E(z) n is included to remove the predicted self-similar 
evolution, where the index n = 4/3, and —1 for the is'(rioou)- 
Tx, i^(riooo)/Jf(0.1r2oo)-7x and Z/x-Tx relations, respectively. 
The scatter expected from statistical uncertainties, <r a tat , is 



(20) 



where 



1 t obs \ 2 , 2 / obs \ 2 

°"i = (fy.i) + a ( cr Tx,iJ > 



(21) 



and o-y^i and cr^f are the errors in T^ s any l^ obs , respectively. 

We estimate the raw scatter, draw, using error- weighted dis- 
tances to the regression line: 



-^^ Wi {log 10 [£(,)T 



(22) 



-alog 10 



TIO 
"X" 



5keV 



log 10 Co 



where w, = cr 2 tat /a 2 and iV° bs is the number of observed NCC 
clusters in the sample. 

Finally, the intrinsic scatter, (Tint, about each observed mean 
relation is estimated as 



2 _ 2 2 
°int — °"raw °"stat • 



(23) 



We examine how well our models reproduce the observed 
scatter about each mean relation, by using the Kolmogorov- 
Smirnov (K-S) test to determine if the residuals for the observed 
and simulated samples are drawn from the same distribution in each 
case. 

Tables [5] and [6] show the probabilities (p- values) for each of 
our 12 models obtained from the \ 2 (equation 1 18t and K-S tests, 
respectively, for the case of the Lx-Tx relation. The models we 
deem to be acceptable are highlighted in bold. It is evident that all 
models with /duty ^ 10~ 3 are ruled out, at least for this particular 
relation, and only three models provide an acceptable match to both 
the mean observed relation and the associated scatter. 

For the sake of brevity, we do not present the corresponding ta- 
bles for the K(riooo)-Tx and A'(riooo)/-R'(0.1r2oo)-Tx relations. 
We simply note that the \ 2 test for the K(riooo)-Tx relation rules 
out all models expect the two with (/rad, /duty) = (1, 10~ 3 ) and 
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Table 6. K-S test probability values for th 
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(1, 10 2 ), whereas the K-S test rules out all models with /duty 
10 ' or / dut y = HT 1 . In the case of the Jsr(riooo)/K(0.1r2oo)- 
Tx relation, both the \ 2 and K-S tests rule out all models with 
/duty = 10 - , tending to favour models that populate the upper- 
right corner of the table. 

We have combined the results of all 6 tests (2 for each of the 
3 scaling relations) using Fisher's method for combining p-values. 
The overall probabilities for our 12 models are summarised in Table 
[7] It is evident that only one model is now acceptable, the model 
with (/ ra d, /duty) = (1, 10~ 2 ). Therefore, we choose our fiducial 
AGN feedback scheme to be the stochastic model JET_R1_D10 -2 . 

We emphasise that the purpose of this section has not been 
to conduct a rigorous statistical analysis, but merely to provide us 
with an indication of the region of our model parameter space that 
is favoured by the observational data. There are several possible 
caveats to our analysis. For example, we are assuming Gaussian 
errors and that each test is strictly independent, both of which may 
not be the case. 



5.3 Metallicity profiles 

Powerful AGN outflows are an obvious candidate for transport- 
ing metal-rich material away from the central regions of ha- 
los. Indeed, there is observational evidence that enriched gas 
is entrained by bubbles inflated by central AGN and removed 
from cluster cores (e.g . iForman et alj|2005l ; iMillion et al.ll2010l : 
Kirkp atrick et alj 1201 lh . Cosmological hydrodynamical simula- 
tions have also demonstrated that A GN are important for the metal 
enrichment of intracluster gas (e.g. ISiiacki et alj|2007fc iMoll et all 
|2007| ; iFabian et alfeOlOl : IWiersma et alfcOllh . 

Recall from Section 1331 that SN feedback alone has a negli- 

Table 7. Combined probability values for all three scaling relations. 
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Figure 11. Emission- weighted metallicity profiles for cluster CI resimu- 
lated with our fiducial stochastic AGN feedback model, but with differ- 
ent metal enrichment schemes (solid coloured lines; see legend for model 
details). The profile obtained from a ran with supernova feedback only is 
shown for comparison, as well as the profiles of observed CC (dashed grey 
lines) and NCC (solid black lines) clusters from the sample of MAT1 1. 

gible impact on the distribution of metals in the ICM, regardless 
of the way in which the available energy is injected into the dif- 
fuse gas. We now investigate whether the inclusion of extra en- 
ergy input from AGN affects ICM metallicity profiles. To do this, 
we have resimulated cluster CI with our fiducial AGN feedback 
model (JET_R1_D10 -2 ), varying the parameter /z.rad in our metal 
enrichment scheme (recall that this parameter sets the radius of the 
spherical region about a galaxy within which metals are injected). 

Figure [TT] compares the resulting emission-weighted metal- 
licity profiles with those of observed clusters from the sample of 
MATH. As before, the best match to the gradient of the profiles 
of observed NCC clusters arises when metals are distributed uni- 
formly throughout the entire halo (/z, r ad — 1)- As /z.rad is de- 
creased, we see the development of a sharp central abundance peak 
that is in conflict with the observational data. Recall that a dis- 
tributed metal enrichment model is justified since the SA model 
underpinning our simulations predicts that almost all metals in the 
ICM are accreted, rather than being produced by BCGs. 

Comparing the metallicity profile predicted by model JET_Z1 
to that obtained from the SN_KiR_Zl run, we can see that addition 
of AGN feedback has indeed displaced some metals from central 
cluster regions, leading to a flatter profile, but the effect is small. 
This is because only a small fraction of the particles in the core are 
heated by AGN in our model, and they receive a sufficiently large 
entropy boost to escape the central regions of clusters, leaving the 
majority of metal-rich material behind in the core. 

5.4 Summary 

We have developed a new AGN feedback model where heat energy 
is injected into intracluster gas in a stochastic manner. Our model is 
physically-motivated and has just two free parameters: / ra d, which 
sets the radius of the spherical region within which energy is in- 
jected, and /duty, which is the fraction of particles in this region 
that are heated per AGN duty cycle. 



We have found that ICM entropy profiles and the Lx-Tsi scal- 
ing relation are fairly insensitive to variations in / ra d, but depend 
strongly on /duty For small values of /duty, the resulting entropy 
profiles are close to those obtained from a run with SN feedback 
only, implying that AGN heating has little effect in this case. This 
is because only a few particles are heated, and they receive large 
entropy boosts, which causes them to rise buoyantly to large dis- 
tances from the cluster centre, leaving the majority of low-entropy 
gas behind. 

As /duty is increased, the probability of a particle being heated 
also increases: hence more particles in the core are heated, they 
are given a smaller entropy injection, and so they do not escape 
central cluster regions. As the heated gas expands, the gas density 
drops, causing the gradient of the resulting entropy profiles to be- 
come shallower. 

We have used three observed scaling relations to identify 
an optimal choice for the two free parameters in our model: 
(/rad, /duty) = (1,10~ 2 ). Setting / d uty = 10~ 2 roughly corre- 
sponds to a jet opening angle of 16° . With these parameter choices, 
our model, named JET_R1_D10~ 2 , can explain the observed en- 
tropy profiles and Ly^-T relations for NCC clusters. 

Using JET_R1_D1(T 2 as our fiducial AGN feedback model, 
we have demonstrated that AGN heating has little impact on the 
distribution of metals in the ICM by comparing to a model with ki- 
netic supernova feedback only. When the metals produced by stars 
in galaxies are distributed uniformly throughout the entire host halo 
(model JET_Z1), the resulting abundance gradients provide a good 
match to those observed in NCC clusters. 

For the remainder of this paper we thus adopt model JET_Z1 
as our fiducial model for star formation, metal production, black 
hole growth and associated stellar and AGN feedback. 



6 COMPARISON WITH OBSERVATIONS 

In this section we conduct a detailed assessment of how well our 
fiducial model (JET.Z1) can reproduce key observed thermal and 
chemical properties of intracluster gas. 

6.1 Thermal properties of the ICM 

Figure[72lcompares the predicted entropy profiles of all 25 clusters 
(C1-C25) in our sample with the profiles of observed systems in 
the same mass range. The profiles of relaxed (disturbed) simulated 
clusters are shown as solid (dotted) red lines. 

First impressions are that our fiducial model generates clus- 
ters whose entropy profiles agree well with those of observed NCC 
systems, both in terms of normalisation and gradient. The central 
entropy is too high in three of our objects, but two of these are 
classified as disturbed systems. To assess our model more quantita- 
tively, we now examine how the entropy profile normalisation and 
shape scale with system temperature. 

Figure[T3]shows the entropy profile normalisation (defined as 
the entropy at riooo) as a function of spectroscopic-like tempera- 
ture. Filled (open) circles correspond to relaxed (disturbed) objects 
(both here and in all subsequent figures). Observational data for 
NCC clusters in the REXCESS are also shown. The parameters 
of the accompanying predicted and observed best-fit relations are 
summarised in Table[8] Note that we only consider relaxed systems 
in our sample when performing the fit and, as before, we adopt the 
BCES orthogonal fitting method. It is evident that the K (riooo)-T B \ 
relation predicted by our model is a good match to the observed 
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Figure 12. Entropy profiles for 25 clusters resimulated with our fiducial 
stochastic AGN feedback model. The profiles of relaxed (disturbed) sys- 
tems are shown by solid (dotted) red lines. Observed profiles of CC (dashed 
grey lines) and NCC (solid black lines) clusters in the REXCESS sample 
(PAP10) are also displayed for comparative purposes. The observed clusters 
span the same mass range as our simulated ones. 



Figure 14. The scaling of entropy profile shape, K (t-lqqq)/ K(0.1r2oa)> 
with temperature predicted by our fiducial stochastic AGN feedback model. 
The spectroscopic-like temperature is computed within r'500 ■ Relaxed (dis- 
turbed) systems in our simulated cluster sample are shown as filled (open) 
circles. The solid red line is the best-fit relation obtained using our relaxed 
clusters only. For comparison, we also display the observed relation for 
NCC clusters from the REXCESS (PAP10; open squares and black line). 
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Figure 13. The scaling of entropy profile normalisation, i^(riooo). with 
temperature predicted by our fiducial stochastic AGN feedback model. We 
compute the spectroscopic-like temperature within rsoo ■ Filled (open) cir- 
cles correspond to relaxed (disturbed) systems in our simulated cluster sam- 
ple, and the solid red line is the best-fit relation considering relaxed sys- 
tems only. The observed relation for NCC clusters in the REXCESS sample 
(PAP 10) is also shown by open squares and a solid black line. 

relation: the slope and normalisation are both within la of the ob- 
served values, and the scatter about the mean relation is compara- 
ble. 

In Figure [14] we display the variation of the ratio 
K (r 1000) /K(0.1r 200) (a measure of the entropy profile shape) 
with temperature. The slope of our predicted relation is consis- 
tent with that of the observed relation and the scatter is iden- 
tical to the observed value; see Table [8] However, the normal- 



isation is slightly lower. There are two reasons for this offset. 
First, one of our relaxed clusters has an anomalously low value of 
K(riooo)/ K(0.1r2oo), which lowers the normalisation of the pre- 
dicted i\~(riooo)/^'(0.1r2oo)-T'si relation. This objects correspond 
to the relaxed system with an excessive central entropy in Figure 
[T2l Second, three of the observed clusters lie considerably above 
any of our simulated clusters on the K(riQQo) / K(0.1r2oa)-T s \ 
plane. This acts to raise the normalisation of the observed relation 
relative to the predicted one. Although classified as NCC systems 
in REXCESS, these objects actually have a low central entropy, 
reminiscent of CC clusters; see Figure [12] Without these outliers, 
there is good overall agreement between the predicted and observed 
^(nooo)/^(0.1r 2 oo)-T s i relations. 

Finally, we contrast our predicted Lx-T s \ scaling relation 
with the REXCESS NCC cluster relation in Figure [B] We pre- 
dict slightly less scatter about the mean relation than observed, but 
we recover the normalisation and slope of the observed relation to 
within la, as summarised in Table [8] We conclude that our fidu- 
cial model yields an Lx-Tsi relation that is a good match to the 
observed relation for NCC clusters. 



6.2 Chemical properties of the ICM 

The predicted emission-weighted metallicity profiles of all 25 clus- 
ters in our sample are displayed in Figure [16] along with Fe abun- 
dance profiles of observed CC and NCC clusters from the sample 
of MATH. To ensure a fair comparison, we only plot the profiles 
of observed clusters that lie in the mass range spanned by our sim- 
ulated sample. Again, solid (dotted) red lines correspond to relaxed 
(disturbed) systems. We remind the reader of the limitations of our 
metallicity model: we assume only prompt enrichment and cannot 
discriminate between ejecta from core collapse and Type la SN. 
Nevertheless, it is evident that the profiles of our simulated clus- 
ters are in reasonable agreement with those of observed NCC clus- 
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Table 8. Best-fit parameters (with lcr errors) for 2 = scaling relations obtained from our full 25-cluster 
simulated sample, and from the REXCESS observations of PAP10. Note that we only consider relaxed clusters 
in our sample when deriving predicted relations. All fits were performed using the BCES orthogonal regression 
method. 
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1013 ± 41 
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0.032 


K(nooo)/K(0.1r2oo)-T s i 
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0.57 ±0.33 
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0.071 
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3.22 ±0.12 


0.098 


Z(0.25n 80 )/Z(0.045ri8o)-T sl 


0.705 ±0.019 


-0.071 ±0.086 


0.055 


0.824 ± 0.067 


-0.11 ± 0.14 


0.029 



Co and a are the best-fitting normalisation and slope of the relations, respectively (see equation [T9l . and (Ti n t is 
the intrinsic scatter about the mean relation (equation [23}. 
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Figure 15. The X-ray luminosity-temperature scaling relation predicted by 
our fiducial stochastic AGN feedback model. X-ray properties are computed 
within T'5oo ■ Filled (open) circles correspond to relaxed (disturbed) systems 
in our simulated cluster sample, and the solid red line is the best-fit relation 
for relaxed objects only. For comparative purposes, we also plot the ob- 
served relation for NCC clusters in the REXCESS sample (PCA09; squares 
and solid black line). 



ters, both in terms of normalisation and slope, with the exception 
of the few observed NCC objects that have a sharp central abun- 
dance peak, which could be CC remnants jLeccardiet"al]|20iCl : 
iRossetti & Molendil l2010). To demonstrate this more rigorously, 
we now investigate how the metallicity profile shape scales with 
temperature. 

The measure of metallicity profile shape we adopt is the ra- 
tio of the metallicity at a radius of 0.25ri8o to that at a radius of 
0.045ri8o. We chose these particular radii since nearly all of the 
clusters in the sample of MATH have a metallicity profile defined 
over this radial range, thereby maximising the number of observed 
clusters we can compare our predictions to. 

Figure [T7] shows the predicted scaling of the ratio 
Z(0.25ri8o)/Z(0.045ri8o) with spectroscopic-like temperature. 
The corresponding observed relation for NCC clusters is also 
shown for comparison, and the parameters of both best-fit relations 
are presented in Table[8] We recover the observed gradient to within 
lcr, but the predicted normalisation is lower than observed, and the 
scatter is larger. However, the observational errors are large and 



Figure 16. Emission-weighted metallicity profiles for 25 clusters resim- 
ulated with our fiducial stochastic AGN feedback model. The profiles of 
relaxed (disturbed) systems are shown by solid (dotted) red lines. For com- 
parison, we also show observed profiles of CC (dashed grey lines) and NCC 
(solid black lines) clusters in the REXCESS sample (PAP10). We only show 
observed clusters with a mass in the same range as our simulated objects. 



there is one observed cluster that lies considerably above all the 
others, which will act to increase the normalisation of the observed 
relation relative to that of ours. 



7 INCLUDING RADIATIVE COOLING: A FIRST 
ATTEMPT 

None of the simulations presented in this paper thus far incorporate 
cooling processes. In this section, we make a first attempt to extend 
our hybrid feedback scheme by allowing gas to cool radiatively. 
Our aim is to formulate a feedback model which can produce both 
CC and NCC clusters, whilst avoiding catastrophic over-cooling 
of gas in central cluster regions. We emphasise that this work is 
exploratory, intended merely to demonstrate that such a model is 
possible with our approach. 

The addition of gas cooling is likely to lead to differences be- 
tween the predictions of kinetic and thermal feedback schemes that 
have not been apparent in our previous non-radiative runs. How- 
ever, it is not our intention here to conduct an exhaustive com- 
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Figure 17. The scaling of emission-weighted metallicity profile shape, 
•Z(0.25ri8o)/Z(0.045ri8o), with temperature predicted by our fiducial 
stochastic AGN feedback model. The filled (open) circles correspond to re- 
laxed (disturbed) systems in our simulated cluster sample. The solid red 
line is the best-fit relation obtained when considering just the relaxed sys- 
tems in our sample. The observed relation for NCC clusters in the sample 
of MATH is also displayed (squares and solid black line). 



parison of different feedback models when cooling processes are 
included; we save this for future work. 

We implement AGN feedback using the 2-parameter stochas- 
tic heating model developed in the Section[5](however, as we shall 
see, the optimal parameter choices change with the addition of 
cooling), and we adopt our fiducial model for SN feedback (gas 
particles neighbouring a galaxy are given a kick in a random di- 
rection with a speed of 600 km s _1 ) and metal enrichment (metals 
produced by stars in galaxies are uniformly distributed throughout 
the entire host halo). 

Metal-dependent radiative cooling is included in our simula- 
tions as follows. For each gas particle, we know its (smoothed) 
metallicity, Z snlt i (see equation [T3l>. and we can compute its tem- 
perature from its entropy, Ai, and density, pi. With this informa- 
ti on we then calculate the coo ling rate using the cooling function 
of ISutherland & Dopitd Jl993l) . and reduce the entropy of the gas 
particle accordingly. 

Figure [T8lcompares the entropy profile of cluster CI obtained 
from runs with our fiducial feedback model (JET_Z1) with and 
without cooling. It is apparent that the heating from SN and AGN 
has not been sufficient to prevent over-cooling in central cluster 
regions: there is a sharp drop in gas temperature at r < 0.3rsoo, 
leading to a steep decline in the entropy profile, and there is an en- 
tropy increase at larger radii due to hotter, lower-density gas flow- 
ing inwards from cluster outskirts to maintain pressure support in 
the core. 

A priori, there is no reason to expect that the amount of SN 
and AGN heating provided by the underlying SA model would be 
sufficient to precisely balance radiative cooling in the simulation. 
This is because L-Galaxies employs a simple cooling recipe based 
on the assumption that haloes have a spherically-symmetric isother- 
mal gas distribution, which is typically not the case in hydrodynam- 
ical simulations, so the predicted cooling rate of gas in haloes will 
be different to that in the simulation. Since the amount of gas that 
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Figure 18. Entropy profiles for cluster CI resimulated with our fiducial 
stochastic AGN feedback model without cooling (solid red line), and with 
metal-dependent radiative cooling (solid blue line). Observed profiles of 
similar-mass CC (dashed grey lines) and NCC (solid black lines) clusters 
in the REXCESS sample (PAP10) are also displayed for comparative pur- 
poses. 



can cool to form stars and accrete onto central black holes governs 
the level of subsequent feedback, such differences in gas cooling 
rates imply that it is unlikely a self-regulating feedback loop would 
be established in the simulation. 

To address this problem, we have developed an ad hoc ex- 
tension of our stochastic AGN feedback model where we inject 
extra energy into cluster cores as a crude representation of addi- 
tional AGN heating that would have arisen from enhanced black 
hole accretion due to more efficient cooling. Such a scheme is justi- 
fiable, provided the extra energy input required to balance radiative 
cooling is a small fraction of that originally available from the SA 
model. 

The details of our model are as follows. At each SA model 
output, we identify all gas particles in the simulation residing in 
the central regions of haloes (r < 0.1r2oo)- At each subsequent 
timestep, we test if any of these particles have cooled below a 
threshold temperature of 3 x 10 4 K; if they have, we raise their 
temperature to some multiple, /temp, of the virial temperature, T200 
(equation|3j, of their host halo at the previous output time. We con- 
tinue in this fashion until the next model output is reached, at which 
point the list of particles contained in halo cores is reset and the 
process is repeated. 

In what follows we keep the radius of energy injection in our 
AGN feedback model fixed at unity (/ ra d = 1)- We then have two 
free parameters: the fraction of particles heated per AGN duty cy- 
cle, /dut y! and /temp, which controls the temperature cold particles 
in cluster cores are heated to. We now explore the effect of vary- 
ing these parameters on the entropy distribution in clusters. All of 
our models are summarised in Table [9] where we have given them 
the label ZCOOL to emphasise that they include metal-dependent 
radiative cooling. 
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Table 9. Stochastic AGN feedback models with metal-dependent radiative cooling and a prescription for addi- 
tional heating of cold gas in cluster cores. In each case, supernova feedback is implemented using a kinetic model 
where particles neighbouring a galaxy are given a kick in a random direction with velocity 600 km s" 1 . Energy 
and metals are both injected within a radius of T200 (/rad = 1 an d /z,rad = 1. respectively). Unless otherwise 
stated, the fraction of particles heated per AGN duty cycle is /duty = 10" 1 , and cold particles in cluster cores 
are heated to a temperature of /temp = 2.5 times the halo virial temperature. 



Model name 


Type 


Energy injection method 


Comments 


ZCOOL_D/ duty 


Stochastic 


Fixed entropy 


/duty = Kr 4 , 10- 3 , 10- 2 , HT 1 


ZCOOL.T/temp 


Stochastic 


Fixed entropy 


/temp = 1,2.5,4.5,20 
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Figure 19. Entropy profiles for cluster CI resimulated with our stochastic 
AGN feedback model, including metal-dependent radiative cooling, for var- 
ious choices of the parameter /duty (coloured lines; see legend for model 
details - note that the exessive cooling for /duty = 10~ 4 has led to reclas- 
sification of the cluster as having high substructure). We also assume that 
any cold gas remaining in cluster cores is heated to 2.5 times the halo virial 
temperature. For comparison, we also show the profiles of observed CC 
(dashed grey lines) and NCC (solid black lines) clusters in the REXCESS 
sample (PAP10). 

7.1 The effect of changing /duty 

Figure [19] illustrates how the entropy profile of cluster CI is af- 
fected by varying /duty, keeping /temp fixed at 2.5. With the excep- 
tion of /duty = 10 _1 , all values of /duty produce entropy profiles 
that exhibit signs of over-cooling. This is because only a small frac- 
tion of core particles are heated by AGN feedback in these cases, so 
gas can cool efficiently in the cluster core, even with the injection of 
additional energy. When we increase /duty to 10~ 1 , so that a larger 
fraction of core particles are heated by AGN, we obtain a CC-like 
entropy profile and it appears that radiative cooling has been bal- 
anced. Therefore, we change our fiducial value of /duty from 10 -2 
to 10 - , which corresponds to a larger jet opening angle of about 
52°. 



7.2 The effect of changing /t C m P 

The effect of varying the parameter /temp on the entropy profile of 
cluster CI is shown in Figure [20] The values of /temp we consider 
are 1, 2.5, 4.5 and 20, fixing /duty = 0.1 in each case. Cluster CI is 



Figure 20. Entropy profiles for cluster CI resimulated with our stochas- 
tic AGN feedback model, including metal-dependent radiative cooling, and 
assuming that cold gas in the central regions of clusters is heated to some 
multiple, /temp, of the halo virial temperature (solid, coloured lines; see 
legend for model details). The other model parameter, /duty, is fixed at 
10 _1 . The profiles of observed CC (dashed grey lines) and NCC (solid 
black lines) clusters in the REXCESS sample (PAP10) are also displayed 
for comparison. 

a relaxed system that has not recently undergone any major mergers 
(half of its mass was in place at z ~ 0.8), so it is a prime candidate 
for developing a CC. It is evident that, as /temp is decreased from 
20 to 1, the entropy profile steepens, becoming progressively more 
like that of a CC cluster. This trend can be explained as follows. 
For small values of /temp, cold particles in cluster cores that have 
received an additional energy input are able to radiate away this 
energy more quickly than when /temp is large, because they have 
not been heated to such a high temperature and thus their cooling 
time is shorter. Therefore, as /temp is decreased, the amount of 
cool, dense gas in cluster cores increases, leading to a lower central 
entropy and a steeper profile. 

As mentioned above, an important issue to address is whether 
our heating model is energetically plausible. To quantify this, we 
define /energy as the ratio of the amount of extra heat energy sup- 
plied to that originally available from the SA model over the course 
of the simulation. We want /energy to be as small as possible at 
z = 0. For /temp = 1, 2.5, 4.5 and 20, we have /energy ~ 0.6, 
0.15, 0.18 and 0.5, respectively, so we discard the models with 
/heat = 1 and 20 on energetic grounds. Small values of /heat 
(/heat = 1) lead to a large extra energy input because heated 
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Figure 21. The X-ray luminosity-temperature scaling relation predicted by 
our fiducial stochastic AGN feedback model with metal-dependent radiative 
cooling and additional heating of cold gas in cluster cores. X-ray properties 
are computed within rsoo- CC (NCC) clusters in our sample are shown 
as blue (red) circles, while filled (open) symbols denote that a cluster is a 
relaxed (disturbed) system. For comparative purposes, we also plot observa- 
tional data for CC (diamonds) and NCC (squares) clusters in the REXCESS 
sample (PCA09). 

gas is able to cool down relatively quickly in core regions, and is 
then heated again, so many extra heating events are required over 
the formation history of a cluster. Conversely, when /heat is large 
(/heat = 20), few extra heating events are required because any 
cold gas is heated to such high temperatures that its cooling time 
becomes very long. However, the large amounts of energy needed 
to heat gas to such high temperatures mean that /energy is again 
large. 

For the remainder of this section, we choose the model with 
(/rad, /duty, /heat) = (1, 10 _1 , 2.5) as our fiducial model since 
this yields a CC-like entropy profile for cluster CI, yet only re- 
quires an extra energy input of ~ 15% of that available from the 
SA model. We have simulated all 25 clusters in our sample with 
this model, and we now examine the predicted thermal properties 
of the ICM. 

7.3 Thermal properties of the ICM 

Figure[2"Tlshows our X-ray luminosity-temperature relation, where 
both luminosity and spectroscopic-like temperature have been 
computed within rsoo- Blue (red) circles represent CC (NCC) clus- 
ters, and filled (open) symbols correspond to relaxed (disturbed) 
systems. We classify objects as CC clusters if they are scattered 
above the mean observed relation for NCC clusters in the REX- 
CESS (PAP10) by more than la; 7 of our 25 objects satisfy this 
criterion. For comparative purposes, we also show observational 
data for CC (diamonds) and NCC (squares) clusters in the REX- 
CESS. 

Our predicted relation is not a perfect match to the observa- 
tional data, in the sense that the slope appears steeper and there are 
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Figure 22. Mean entropy profiles for clusters. The upper, red lines corre- 
spond to NCC systems, and the lower, blue lines to CC systems. The solid 
lines are model clusters resimulated with our fiducial stochastic AGN feed- 
back model, plus metal-dependent radiative cooling and additional heating 
of cold gas in cluster cores. The dashed lines are observed profiles of clus- 
ters in the REXCESS sample (PAP10) that straddle the same mass range. 

no low-temperature systems with a high luminosity, although this 
could be a selection effect. However, the salient point is that we 
are able to generate both CC and NCC systems with a single feed- 
back model, a feat that is notoriously difficult with self-consistent 
hydrodynamical simulations. This encouraging result warrants fur- 
ther development of our model in future work. 

The mean entropy and spectroscopic-like temperature profiles 
of all 25 clusters in our sample are displayed in Figures [22] and [23] 
respectively. The profiles of CC (NCC) clusters are shown by blue 
(red) lines, with solid (dashed) lines corresponding to simulated 
(observed) systems. Note that the dispersion of individual clusters 
about these mean relations is quite large, especially for the temper- 
ature profiles, so that it would not be possible to look at a paricular 
profile and classify it with certainty as either NCC or CC, according 
to our definition above. 

There are a number of diferences between the simulated and 
observed profiles. Firstly, both the NCC and CC simulated temper- 
ature profiles are too low in the cluster cores. The entropy profile of 
simulated CC clusters has the correct slope, but too high a normal- 
isation below 0.3 r-500- Finally, the simulated NCC entropy profile 
shows no sign of flattening at the smallest radii. These features all 
suggest that our heating model is far from perfect and that per- 
haps we should target additional heating at not just the coldest gas. 
Nonetheless, it is pleasing that this first attempt should at least lead 
to a distinct separation in the mean profiles of the two classes of 
cluster. 

Finally, we assess the energy requirements of our model. In 
Figure[24] we show the ratio /energy as a function of redshift, aver- 
aged over all 25 clusters. The solid (dotted) lines show the differen- 
tial (cumulative) evolution, and again we have divided our sample 
into CC (blue lines) and NCC (red lines) systems. The first point 
to note is that, at z = 0, the average total extra energy input is 
~ 15% of that available from the SA model (this is actually true 
for all but one of our objects, which has /energy ~ 0.25), which is 
reasonable. CC systems require a larger total energy input, which 
is to be expected since the gas cooling time in the central regions 
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Figure 23. Mean spectroscopic-like temperature profiles for clusters. The 
upper, red lines correspond to NCC systems, and the lower, blue lines to CC 
systems. The solid lines are model clusters resimulated with our fiducial 
stochastic AGN feedback model, plus metal-dependent radiative cooling 
and additional heating of cold gas in cluster cores. The dashed lines are 
observed profiles of clusters in the REXCESS sample (PAP10) that straddle 
the same mass range. 
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Figure 24. Evolution of the ratio of the additional energy input required 
to offset radiative cooling in cluster cores to that available from SN and 
AGN feedback, averaged over all 25 clusters in our sample. We have split 
our sample into CC (blue) and NCC (red) clusters. The solid (dotted) lines 
show the differential (cumulative) evolution. 



of such objects is shorter than in NCC objects, but the difference 
is small. Interestingly, for redshifts z < 3, there is little to distin- 
guish between the two; the main difference occurs at high redshift, 
z ~ 3 — 8, where, on average, much more additional energy is in- 
jected into the cores of CC clusters than NCC clusters. This may be 
an indication of an earlier assembly history for CC clusters. 



8 CONCLUSIONS 

In this paper we have used numerical simulations to investigate 
how star formation, black hole accretion and the associated feed- 
back from SNe and AGN heat and enrich intracluster gas. Our 
primary objective was to assess how different implementations of 
these feedback processes affect the thermal and chemical proper- 
ties of the ICM, using a selection of data from X-ray observational 
studies to constrain our models. 

We have resimulated a sample of 25 massive galaxy clusters 
extracted from the Millennium Simulation. In these simulations, 
the energy and metal input into the ICM by SNe and AGN is cal- 
culated from a SA model of galaxy formation, using the hybrid 
scheme of SHT09. This guarantees that feedback originates from a 
realistic galaxy population, whereas fully self-consistent hydrody- 
namical simulations often predict excessive star formation on clus- 
ter scales. 

Our main achievement has been to develop a new model for 
AGN feedback that is both physically motivated and capable of ex- 
plaining several fundamental observational properties of clusters. 
All of the other, more commonplace, models we have tested fail on 
one or both of these points. Our new model is based on stochas- 
tic, anisotropic heating of the ICM, which is motivated by observa- 
tional evidence that AGN heating is likely to be directional, rather 
than isotropic. 

Our conclusions are as follows: 

(i) Energy input from SN-driven galactic winds has no effect on 
the entropy and metallicity structure of the ICM, regardless of the 
method used to inject energy and metals into the intracluster gas. 

(ii) Simple thermal AGN feedback models all heat the gas ex- 
cessively in the central regions of clusters, generating flat ICM en- 
tropy profiles that disagree with the observational data. Differences 
between our various models are negligibly small, even though the 
number of particles heated by AGN can vary enormously between 
models. 

(iii) Kinetic AGN feedback models can reproduce the observed 
entropy profiles of NCC clusters, but only if the wind speed is very 
high, u w i n d = 20 000 km s , which is possibly unreasonable on 
physical grounds. The success of this model is due to the fact that 
only a small number of particles near the centre of the halo are 
kicked, and the momentum boost they receive is sufficient to trans- 
port them to cluster outskirts, leaving low-entropy gas behind in the 
core. 

(iv) There are two free parameters in our stochastic heating 
model: / ra( j, which governs the radius of the region (in units of the 
virial radius) about a galaxy in which energy is injected, and the 
fraction, /duty, of neighbouring gas particles that are heated per 
AGN duty cycle. The parameter /duty can be linked to the open- 
ing angle of AGN jets, assuming jets can be simply modelled by 
biconical outflows. Using the observed scaling of X-ray luminos- 
ity, entropy profile normalisation and shape with temperature as 
constraints, we identified (/ ra d> /duty) = (1, 10 -2 ) as an optimal 
choice for these parameters. The choice /duty = 10~ 2 corresponds 
to a jet opening angle of roughly 16°. This is our fiducial model for 
AGN feedback. 

(v) Our fiducial stochastic heating model is able to explain both 
the thermal and chemical properties of intracluster gas, at least for 
NCC systems. We obtain a good match to several key pieces of 
observational data: the normalisation and shape of entropy profiles, 
the X-ray luminosity-temperature scaling relation, and the shape of 
metallicity profiles. The model is successful for the same reason as 



22 



C. J. Short et al. 



the kinetic AGN feedback model with u w i n( j = 20 000 km s -1 , but 
has the advantage of being physically motivated. 

(vi) Reproducing the observed abundance gradient in NCC clus- 
ters requires that metals ejected from galaxies are distributed 
throughout the entire halo. Injecting metals in a concentrated fash- 
ion leads to a sharp central peak in the metal distribution that is not 
observed. A metal enrichment model where metals are distributed 
throughout the halo is consistent with the SA model underlying our 
simulations, which predicts that over 95% of the metals in diffuse 
gas are accreted, rather than being produced by the central galaxy 
of the halo. 

(vii) AGN heating causes a flattening of ICM metallicity pro- 
files, but the effect is small in our fiducial model. This is because 
only a small fraction of the particles in the core are heated by AGN, 
and they receive a sufficiently large entropy boost to escape the cen- 
tral regions of clusters, leaving the majority of metal-rich material 
behind in the core. 

(viii) With the addition of metal-dependent radiative cooling, 
our stochastic AGN feedback model is capable of producing CC 
and NCC systems, avoiding catastrophic over-cooling, but only if 
we assume that additional energy is injected into cold gas in clus- 
ter cores to offset radiative losses. The justification for this simple 
model is that we expect a mismatch between gas cooling rates in 
the SA model and hydrodynamical simulations. The amount of ex- 
tra energy typically required is ~ 15% of that available from SN 
and AGN feedback over the formation history of a cluster. 

(ix) As Figures [22] and [23] illustrate, our cooling model is far 
from a perfect match to the observed entropy and temperature pro- 
files in cluster cores. Nonetheless, they do show a distinction be- 
tween CC and NCC clusters. 

To keep the model for the ICM as simple as possible, we have 
neglected a number of physical processes in this paper. While we 
do not expect any of these to make a major contribution outside the 
core of the cluster, it is possible that their cumulative effect could be 
important. This should be investigated further in future extensions 
of this work: 

• Magnetic fields and cosmic rays will help to provide ex- 
tra pressure support in the ICM. They can be generated both by 
mergers and by AGN activity. However, observations suggest that 
neither makes a domina nt contribution to clusters, excpet per- 
haps in the core regions. iBrunettil fcOllh summarises the current 
state of play for cosmic rays: gamma ray observation s from Fermi 
dAckermann & et alj|2010l: Ijeltema & Profumol F20 1 ) limit the en- 
ergy density of comsic rays to less than a few hundredths of that of 
the thermal energy of the ICM. We note that the qualitative effect 
on the cluster gas den sity profile of the inclusion of cosmic rays 
differ s between AMR dVazza et al]|2012h and SPH Jjubelgas etail 
2008) simulations, but the effect on the density and temperature 
profiles of clusters is minor in each case. 

Magnetic field strengths in the cores of cluste r s range from a few 
to a few tens of fj,G (e.g. iBonafede et"aHl20ld : IVacca et alj|2012l 
and references therein). Measurements for the cluster as a whole 
are hard to make and generally involv ed a degree of model ling. 
The observations have been reviewed by Bonafede ~et~ai] l l2010l) and 
indicate typi cal values o f J-2uG. This agrees with theoretical es- 
timates from lKunz et al. 

I {III 

]]). At this level, the magnetic energy 
density will be only a minor contributor (of order a percent) to the 
total energy density of the ICM. 

• Conduction has been investigated by many authors, princi- 
pally as a way of overcoming the over-cooling problem in cluster 
cores (e.g. iGuo & Ohl2009l : liparrish et alj201fj| : lRuszkowski & Ohl 



l2010l : iRuszkowski & Ohl l2010l : IParrish et alj |2012|) . Direct evi- 
dence for conduction is however, by its very nature, almost impossi- 
ble to achieve. All we have are upper limits based on the existence 
of large temperature gradients surrounding clumps of hot or cold 
gas (e.g. iGu et alJ2009llde Plaa et alj20iq|McDonald et alj2012l: 
iRussell et all2012D . It is possible that these clumps are surrounded 
by magnetic sheaths that limit conduction across the interface: for 
the purposes of this paper, what is important is the degree by which 
magnetic fields would suppress large-scale conduction between the 
core and the cluster outskirts, and this will depend upon the rela- 
tive importance of ordering by convective motions and the stirring 
by galaxies and infalling substructure. 

Using the observ ed density and temperature profiles of 
Vikhl ininetalji 2006), and assuming conduction at the Spitzer rate, 
then it is possible to estimate the maximum rate at which gas can 
be heated or cooled. This can be quite large in the cluster core, 
but is of order (lOGyr) for r > O.lrsoo (for a 3keV cluster). 
Given that some suppression below the Spitzer rate is likely, then 
conductive heat transport at these radii will be minor, but perhaps 
not completely negligible. We note that the effect would be to heat 
gas within the core and cool gas at larger radii, thus flattening the 
entropy gradient even further and reinforcing the arguments in this 
paper. 

We have demonstrated that our fiducial stochastic heating 
model can explain several important observational properties of 
massive clusters, at least for those systems without an X-ray bright 
CC. With the inclusion of metal-dependent radiative cooling and 
a simple prescription for additional heating of gas in the central 
regions of clusters, we have taken our first steps to being able to 
produce both CC and NCC systems with a single model. We are 
currently undertaking a hydrodynamical simulation of the full Mil- 
lennium volume (500/i -3 Mpc 3 ) with these models. The aim is to 
produce a large, publicly-available sample of galaxy groups and 
clusters whose properties are consistent with the available X-ray 
data. An example of an important application of such a sample 
would be modellin g the selection functions of X-ray surveys (e.g. 
ISahlen et al]|2009h . This is essential to exploit the full power of 
clusters as cosmological probes of the expansion history of the Uni- 
verse. 
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